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Abstract. The angular motion of a few-body system is described with global 
vectors which depend on the positions of the particles. The previous study 
using a single global vector is extended to make it possible to describe both 
natural and unnatural parity states. Numerical examples include three- and 
four-nucleon systems interacting via nucleon-nucleon potentials of AV8 type 
and a 3a system with a nonlocal era potential. The results using the explicitly 
correlated Gaussian basis with the global vectors are shown to be in good 
agreement with those of other methods. A unique role of the unnatural parity 
component, caused by the tensor force, is clarified in the state of 4 He. Two- 
particle correlation function is calculated in the coordinate and momentum 
spaces to show different characteristics of the interactions employed. 



1 Introduction 

Realistic nucleon-nucleon potentials contain strong tensor components. The ten- 
sor force induces the coupling of different orbital and spin angular momenta in 
the wave function of a few-nucleon system. Because of this property a variational 
solution of the system faces complexities in describing the orbital motion. There 
are two widely used methods to expand the orbital part of the wave function. One 
is a hyperspherical harmonics (HH) expansion [1, 2], and the other is a successive 
coupling of partial waves for the relative motion of the particles [3, 4, 5, 6, 7]. A 
converged solution is attained by increasing the orbital angular momenta. 

Another representation for the orbital motion, proposed in refs. [8, 9], is 
to use a global vector (GV) which is defined by a linear combination of the 
relative coordinates. The coefficients of the combination determine the vector 
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responsible for the rotation of the system. The spherical part of the orbital motion 
is expressed in terms of an explicitly correlated Gaussian (CG). The efficiency of 
the GV representation (GVR) has been tested in many examples [10, 11, 12, 13] 
including not only nuclear but also atomic and hypernuclear physics problems. Its 
accuracy was confirmed in ref. [11] by comparing with the partial- wave expansion 
(PWE) for 3 H and a-particle interacting with central and tensor forces. In these 
studies only a single GV is used, and its application is thus restricted to natural 
parity states. 

The purpose of this paper is to extend the GVR to describe the orbital motion 
with unnatural parity and to test the power of its representation. The unnatural 
parity state is usually unfavorable to gain the energies of low-lying levels, but 
it is needed to reach a precise energy particularly when the noncentral forces 
are taken into account. For example, the ground state of 4 He can be specified by 
three channels with (L, 5)=(0,0), (1,1) and (2,2), where L and S are respectively 
the total orbital and spin angular momenta. Here the channel with L=l and 
positive parity is unnatural, and its mixing probability is of order only 0.37 % [14]. 
However, the L=l unnatural component gives rise to a contribution of more than 
1 MeV to the ground state energy. 

It is not a surprise that the interaction between composite particles like the 
nucleon becomes nonlocal [15], and indeed there is now available a new class of 
nonlocal nucleon-nucleon interactions, e.g., N 3 LO [16] and low-fc [17] potentials. 
The nonlocal operator also appears when one needs to eliminate some redun- 
dant or Pauli-forbidden states from the relative motion between composite par- 
ticles [18]. Thus we take a certain class of nonlocal operators into consideration 
in the present paper. 

A simple method of incorporating the unnatural parity state is to introduce 
two GVs. This was performed in the case of central forces and successfully applied 
to a search for excited states of a positronium molecule [19]. Here we develop a 
method of calculating the matrix elements of not only noncentral forces but also 
other local operators of physical interest such as a multipole density and some 
nonlocal operators. 

The paper is organized as follows. The CG basis with two GVs is defined in 
Sect. 2 and its generating function is introduced. A nice property of the CG basis 
is that its Fourier transform is again the CG in momentum variables. A basic 
method of calculating matrix elements is explained in Sect. 3. The matrix element 
of a nonlocal operator between the generating functions is given in Sect. 4. In 
Sect. 5 we present results calculated for three- and four- nucleon systems using the 
two GVs and compare them to other calculations. As an example of the nonlocal 
potential, we consider a 3a system interacting via the nonlocal a-a potential 
derived from the resonating group method (RGM) [20]. A summary is drawn in 
Sect. 6. Formulas used in the present calculation are collected in Appendices. We 
give the CG basis in the momentum space in Appendix A, present the matrix 
elements of physically important operators in Appendix B, discuss the integral 
transform of the CG basis in Appendix C, and derive the matrix elements for 
some nonlocal operators in Appendix D. 
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2 Correlated Gaussian basis with double global vectors 

The basis function in LS coupling scheme takes a form 

where the square bracket [. . .] stands for the angular momentum coupling. The 
spin and isospin parts are expanded using the basis of successive coupling, e.g., 

Ml = \[---[[[\\\s^\s^Y--]sM s ). (2.2) 

Here the set of intermediate spins (512, 5*123; • • •) takes all possible values for a 
given S. The orbital part V^m"^ ^ s expanded in terms of the CG basis which is 
explained in detail below. The symmetry of identical particles has to be imposed 
on the basis function. 

2. 1 Generating function for correlated Gaussian basis 

Assuming that the system consists of ./V particles, we denote the single-particle 
coordinate of particle i by r; L . It is convenient to introduce relative and center 
of mass coordinates, denoted Xj, to describe intrinsic excitations of the system. 
Both coordinates are related to each other by a linear transformation 

N N 
3=1 3=1 

where x^ is chosen to be the center of mass coordinate. Let x denote an 
(A — l)-dimensional column vector, excluding x^, whose ith element is a usual 
3-dimensional vector aij. A choice of x often employed is the Jacobi set of coor- 
dinates, but may be any other set of relative coordinates. 

A natural parity basis with the orbital angular momentum LM and parity 
(— 1) L is described as 

exp ^-^xAx^j y L M(u\x), (2.4) 

where yLM^uix) is a solid spherical harmonics 

y L M{uix) = \u x x\ L Y LM {u x x) (2.5) 

which is specified by a single GV, U\ x = ^2d=\ u ii x i- The symbol ~ denotes the 
transpose of a matrix. In Eq. (2.4) A is an (TV— 1) X (TV— 1) positive-definite, 
symmetric matrix, and xAx is a short-hand notation for ^2dj~=\ ^ij x i ' x j- The 
matrix A and the (A^— l)-dimensional column vector u\ are parameters to charac- 
terize the "shape" of the basis function. The function (2.4) is a natural extension 
of exp(— \ar 2 )yi m (r) for a single variable case to many-particle system. 

Another extension commonly used is a successive coupling of the partial 
waves corresponding to the respective coordinates, namely, exp(— ^a\xf— ^a2x\— 
■■■) [■■■ [[y h (xi)y i2 {x 2 )}L 12 ye3 ( x 3)]l 123 ■■ -]lm- The basis function (2.4) was 
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compared to that of PWE [11] and the GVR is found to be an excellent al- 
ternative to describe the rotational motion. The angular part was extended in 
refs. [10, 11] to \u\x\ 2K yiMi^ix) , which was found to be advantageous to cope 
with short-ranged repulsive forces. We use the form (2.5) in this paper, however, 
since it greatly simplifies the calculation of matrix elements. 

An unnatural parity basis with the angular momentum LM and parity 
(— 1) L+1 is described using two GVs as 

exp f-^xAx) [yL(y^x)y 1 (u^x)] LM , (2.6) 

where u 2 is an (./V— l)-dimensional column vector which defines the second GV. 
Both forms of Eqs. (2.4) and (2.6) are unified as 

F(l 1 l 2 )lm(ui,u 2 ,A,x) = exp f~xAx\ [yL 1 (uix)y L2 (uix)] LM - (2.7) 

That is, L\ = L, L 2 = for the natural parity case and L\ = L, L 2 = 1 for the 
unnatural parity case. We stress that some important symmetry properties such 
as translation invariance and rotation invariance are already built in the basis 
(2.7). We show in Appendix A that the Fourier transform of Eq. (2.7) is again a 
combination of CG in momentum space. 

The basis function introduced above provides us with a trial function for all 
states but L = and a negative parity. The angular part of the basis for such a 
particular case must contain at least three vectors like Xi ■ (xj x x^) [19]. 

We present calculation formulas of matrix elements for the basis (2.7) with 
arbitrary L\ and L 2 values for a given L. Formulas for some simple operators are 
already given in refs. [8, 9, 10, 12] for natural parity and in ref. [19] for unnat- 
ural parity. In the appendices we give formulas for various operators including 
noncentral potentials as well as those nonlocal kernels which appear in the RGM 
formulation for nuclear cluster models [20]. 

Our method is based on the use of the generating function g for the CG: 

g(s; A, x) = exp ^ — ^xAx + 'sx^j , (2.8) 

where s is an (N — l)-dimensional column vector whose ith element is a 3- 
dimensional vector S{. By expressing with 3-dimensional unit vectors e\ and 
e 2 as Si = \ie\Ui i +\ 2 e 2 u 2i , the basis function (2.7) is generated as follows: 

F( Li l 2 )lm(ui,U2,A,x) = ^yjr^f / / de 1 de 2 [YL 1 (ei)YL 2 (e 2 )]LM 

QLi+L 2 

X dX^dX^ 2 

where 



g(Xie 1 u 1 + X 2 e 2 u 2 ;A,x) , (2.9) 

Ai=A2=0 



B L -Q^. (2 ,0) 

When g(Xie\Ui + X 2 e 2 u 2 ; A, x) is expanded in powers of Ai, only the term of 
degree A^ 1 contributes in Eq. (2.9), and this term contains the Li-th degree 
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of e\ because Ai and e.\ always appear simultaneously. In order for the term 
to contribute to the integration over e±, these L\ vectors of e\ must couple to 
the angular momentum L±, that is they are uniquely coupled to the maximum 
possible angular momentum. The same applies to A2 as well. 

2.2 Coordinate transformation 

Suppose that the coordinate set x is transformed to a new set of relative coordi- 
nates y through y = Vx with an (N— 1) x (N— 1) matrix V. The transformation 
of this kind is needed, e.g., when the permutation symmetry of the constituent 
particles is imposed on the basis function or different coordinate sets are used to 
describe particular correlated motion [4]. The CG basis function in the coordi- 
nate set y is rewritten as 

F( Ll L 2 )LM(ui,u 2 ,A,y) = F( LlL2 ) LM (u Pl ,u V2 ,A p ,x) (2.11) 

with 

A V = VAV, u Vl =Vui, u V2 = Vu 2 . (2.12) 

The coordinate transformation ends up redefining A, u± and u 2 as above, and we 
do not need to introduce different coordinate sets. This property that the func- 
tional form of the CG remains unchanged under the coordinate transformation 
enables one to unify the method of calculating the matrix element. 

3 Calculation of matrix elements 

We introduce the following abbreviated notation for a matrix element 

(F'\0\F) = {F {L3L4)LIMI (u 3 ,u 4 ,A',x)\0\F {LlL2)LM ( Ul ,u 2 ,A,x)), (3.1) 

where O stands for an operator of interest. Using Eq. (2.9) in Eq. (3.1) enables 
one to relate the matrix element to that between the generating functions: 

(F'\0\F) = 

x 

with 

s = Xieiui + X 2 e 2 u 2 , s' = X 3 e s u s + XAe A U4. (3.3) 

The calculation of the matrix element consists of three stages: 

1. Calculate the matrix element, 

M = (g{\3e 3 u 3 +X4e i U4;A',x)\0\g(X 1 e 1 u 1 + X 2 e 2 u 2 ;A,x)), 
between the generating functions. 

2. Expand M in powers of Aj and keep only those terms of degree Lj for each 
i. The terms should contain the unit vector ej of degree Lj as well. One 
may drop any term with A?e»-ei = A? etc. 



IItt/ de ) ([yLs(e3)Y u (e 4 )} LIM/ y[Y Ll (e 1 )Y L2 (e 2 ) 



\LM 



\i=l 



n 

\i=i 



OX,, 



(g(s';A',x)\0\g(s;A,x)) 



Ai=A2=A3=A4=0 



(3.2) 
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3. Perform the differentiation and integration prescribed in Eq. (3.2). 

Formulas for most important matrix elements M are tabulated in ref. [9]. The 
angle integration in stage 3 can be performed using a formula 



11 dei([Y L3 (e 3 )Y L4 (e 4 )WO*[^i(ei)^ 2 (e 2 )]LM 
i=i 

x[[y Ll (e 1 )y L2 (e 2 ) ] L [y L3 (e3)yL 4 (e4)]L']^ 
= ( _ 1) L-LWL.^|+i (iM(t/1 | £ » Jf / )i (3 .4) 

Appendix B collects formulas for the matrix elements of some local operators. 

4 Basic matrix elements for nonlocal operators 

Let us consider the matrix element of a nonlocal operator acting between par- 
ticles k and /. The relative distance vector, r^ — ri, can be expressed as a linear 
combination of Xj (excluding the center of mass coordinate x^) as 

N-l 

Tk ~ Ti = ^ w i x i = WX, (4.1) 
i=l 

where the (N — l)-dimensional column vector w is given by 

Wi = {U- l ) kl - (17-%. (4.2) 
See Eq. (2.3). The nonlocal operator may in general be expressed as 

Vnl = J J drdr'V(r',r)\S(wx-r')){S(wx-r)\, (4.3) 

where V(r', r) determines the form factor of the nonlocal operator. 

The operator of type (4.3) was applied to the study of a-cluster condensation 
in 12 C and 16 [21]. Evaluating the matrix element of the nonlocal operator 
requires deliberation because the coordinate wx has to be singled out from the 
set of (N — l) coordinates Xj, and because the rest of other (iV — 2) coordinates 
must be properly defined on the condition that they are all independent of wx. 
Though the choice of these (N—2) coordinates is not unique, the matrix element 
should not depend on its choice. This problem is discussed in Appendix C. 

As derived in Appendix D.l, the basic matrix element reads 

(g(s>; A', x) \ Vnl | g(s; A, x)) = (^^) * exp Q(5 + 2) J~\ 8 + a')) 

x j j dr'drV(r', r) exp ^— 7^ a ' r ' 2 ~ ^ ar2 ~ fi r> ' r 
x exp ^[(C - rf)s' - rfs] ■ r' + [((-rf)s - rjs f ] ■ r^j , (4.4) 
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where 

B = A + A', c=(wB- 1 W y 1 , J' 1 = B- 1 - cB' 1 ww B' 1 , (4.5) 

and 

a = (A(, a=(A'(, V = J^M, r/ = J^A'C, 

a = a-(Ari, a' = a' -(AW, P = ~C(Arf + A'rj). (4.6) 

Here £ is an (TV— l)-dimensional column vector that is determined uniquely from 
w defined in Eq. (4.1). See Appendix C.2. 

We perform the integration in Eq. (4.4) for a given V(r',r) and follow the 
procedure described in Sect. 3. Some examples are collected in Appendix D. 

5 Numerical results 

5.1 Spectroscopic properties of few-nucleon systems 

A few-nucleon system interacting via a realistic two-body force offers good ex- 
amples of testing the CG basis function with the GVR. The iV-nucleon system 
is specified by the following Hamiltonian 

N N 
i=l j>i=l 

Three-body forces are neglected though their inclusion does not cause any prob- 
lem. The operator Tj is a kinetic energy and the center of mass kinetic energy 
T cm is subtracted in the Hamiltonian. A realistic nucleon-nucleon force is char- 
acterized by a short-ranged repulsion and a long-ranged tensor force. To clarify 
the role of these properties, we employ three types of potential models as Vij. 
The first is the Minnesota (MN) potential [22] which contains rather mild short- 
ranged repulsion and renormalizes the effect of the tensor force into its central 
term. No spin-orbit component of the MN potential is included. The second is 
the AV8' potential [23] which is obtained from the Argonne V18 potential with a 
suitable renormalization procedure. The AV8' potential consists of central, tensor 
and spin-orbit terms. A few-body calculation with this potential shows slow con- 
vergence because of its strong short-ranged repulsion, and in addition the radial 
form factor of AV8' makes the evaluation of matrix elements time-consuming in 
the present approach. The third is the G3RS potential [24] whose radial form is 
given as a combination of Gaussians with three ranges, which makes the numer- 
ical calculation much faster. Other features of the G3RS potential compared to 
AV8' are that the central force is deep, the tensor force is weak and the repulsion 
at the origin is mild. The original G3RS potential contains L 2 and quadratic L S 
terms in even partial waves. The contribution of these terms is small [25], and 
they are ignored in the present calculation. 
The form of v is written as 



vi2 = V c (r) + VbouiMPiA + V t (r)S 12 + V h (r)L • S, 



(5.2) 
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where S12 and LS stand for the tensor and spin-orbit operators and Pi- K = l for 
proton and for neutron. The neutron projection operator Pi u is defined in a 
similar way. The inputs used in this subsection are jra^=^\. 47106 MeVfm 2 
and e 2 =1.440 MeVfm. Here tun is the nucleon mass, and the mass difference 
between proton and neutron is ignored. No isospin mixture is taken into account. 
The u parameter of the MN potential is set to u=l. 

The wave function of the system is expressed as a combination of different 
channel components specified with L and S (see Eq. (2.1)): 



denotes the probability of finding the system in the channel (L,S). Possible 
channels for the ground states of N < 4 systems are (L, S)=(0, 1), (2,1) for 
deuteron, (0,1/2), (2,3/2), (1,1/2), (1,3/2) for 3 H and 3 He, and (0,0), (2,2), 
(1,1) for a-particle. Among these channels, (1,1/2) and (1,3/2) for 3 H ( 3 He) 
and (1,1) for a-particle are unnatural parity. Table 1 compares the deuteron 
properties calculated with the three potentials. Comparing the result between 
AV8' and G3RS, we see that the AV8' potential is stronger in the tensor and spin- 
orbit components but weaker in the central component than the G3RS potential, 
as already noted above. 

The orbital part of the full wave function is expanded in terms of the CG 
(2.7). The CG for the natural parity state has L\=L and ^2=0, and contains the 
parameters A and u\, while the one for the unnatural parity state has L\=L and 
L2=l, an d contains A, u\ and U2- These parameters are selected by the stochastic 
variational method (SVM) [8, 9] as follows: First, we randomly choose a chan- 
nel specified by Lj, L, S from among those (L,S) channels which are included 
in the calculation. For a given S value, a set of intermediate spins is randomly 
chosen. For example, for three-nucleon system with 5=1/2, |[[||]o|]i/2M,g) or 

Table 1. Comparison of deuteron properties calculated with the different potentials. The energy 
and point nucleon root mean square (rms) radius are given in MeV and fm, respectively. The 
probability P(L, S) is given in %. 




(5.3) 



LS 




(5.4) 



Potential 



MN 



G3RS 



AV8' 



E 



-2.202 
10.487 
-12.689 



-2.277 
16.478 
-7.294 
-11.460 



-2.242 
19.881 
-4.458 
-16.641 
-1.024 
1.961 
94.23 
5.77 



(T) 

(Vc) 

(Vt) 
(V b ) 



P(0,1) 

mi) 



1.952 
100 



1.979 
95.22 
4.78 
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| [[55)1 5)1/2 Ms) 1S randomly chosen. A set of intermediate isospins is chosen ran- 
domly as well. After choosing these discrete sets, we determine the nonlinear 
parameters A, m following the SVM procedure. The elements of A are chosen 
randomly from a physically important interval. The diagonal and off-diagonal 
elements of A describe the correlated motion of the system. The parameter m 
is also chosen randomly under the condition U{Ui=\. This parameter serves to 
express the partial waves needed to represent the rotational motion but its role 
is usually not as important as that played by A. Therefore we make more efforts 
to select a suitable A. 

We decrease the energy first by increasing the number of basis functions one 
by one up to a certain dimension. Keeping the number of basis functions fixed, 
we then switch to a refinement process [8, 9] in which each basis function already 
selected is tested against other randomly chosen candidates. We repeat these two 
optimization procedures until a fair convergence is attained. When the number of 
basis functions becomes large enough to get a converged solution, it is convenient 
to rearrange the selected basis set according to importance. By the importance 
we mean the following. Suppose that the number of the basis functions is K. The 
first basis $1 is the one that gives the lowest energy among the K functions. The 
second basis $2 is the one that gives the lowest energy together with $1 among 
the K— 1 functions excluding $1. This ordering is continued until the last basis 
is determined. After rearranging the basis functions in this way, we often 
have those $j which play a very minor role in lowering energy. This is the case 
particularly when i is close to K. These inactive bases may be discarded from 
the basis set to save the basis dimension. After this contraction of the basis set, 
we may again enlarge the basis dimension to search for better basis functions. 

Table 2 lists the energies of three- and four-nucleon systems and the con- 
tributions of the respective terms of the Hamiltonian together with the nucleon 
rms radii. The GVR performance is tested by comparing to other calculations, 
particularly PWE calculations. The basis function in PWE is also Gaussian but 
no explicit correlated terms between the different coordinates are included, so 
that the matrix A is always chosen to be diagonal. A correlated motion is, how- 
ever, accounted for by expanding the trial wave function in different coordinate 
sets. For example, both coordinates of K- and H-types are employed to obtain 
the solution for a-particle. Thus this PWE calculation is the same as the SVM of 
ref. [14], and its accuracy is well tested. A noteworthy difference between GVR 
and PWE is that the latter expresses the angular part of the wave function by 
successively coupling the partial waves ii to the resultant L. In the present cal- 
culation ti is restricted to i± < 2. The total energy of 3 H calculated with AV8' 
and its decomposition to each term agrees with the results of PWE and Fad- 
deev [26]. The agreement between the GVR and PWE calculations is also very 
good for the properties of the a-particle. The four-nucleon system is solved by 
the different methods [14] using the AV8' potential with the Coulomb force being 
neglected. If the Coulomb contribution of the present calculation is omitted, the 
GVR energy becomes —25.85 MeV, which agrees with the lowest energy of the 
benchmark calculations within 70keV. These results confirm that the CG basis 
with GVR is versatile enough to provide us with such accurate solutions that 
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Table 2. The ground state properties of three- and four-nucleon systems. The energy and rms 
radius are given in MeV and fm. The probability P(L, S) is given in %. The basis dimensions 
for 4 He are 400 (G3RS) and 600 (AV8') in GVR, and reduced to half in PWE, respectively. 



Potential 


MN 


G3RS 




AV8' 




Method 


GVR 


GVR 


PWE 


GVR 


PWE 


Ref. [26] 


3 H(I + ) 


E 


-8.38 


-7.73 


-7.72 


-7.76 


-7.76 


-7.767 


(T) 


27.21 


40.24 


40.22 


47.59 


47.57 


47.615 


(Vc) 


-35.59 


-26.80 


-26.79 


-22.50 


-22.49 


-22.512 


(V) 




-21.13 


-21.13 


-30.85 


-30.84 


-30.867 


(V b ) 




-0.03 


-0.03 


-2.00 


-2.00 


-2.003 


vV> 


1.71 


1.79 


1.79 


1.75 


1.75 




P(0, 1/2) 


100 


92.95 


92.94 


91.38 


91.37 


91.35 


P(2,3/2) 




7.01 


7.02 


8.55 


8.57 


8.58 


P(l,l/2) 
P(l,3/2) 




0.03 
0.02 


0.03 
0.02 


0.04 
0.02 


0.04 
0.02 


}0.07 


3 He(i + ) 


E 


-7.71 


-7.08 


-7.08 


-7.10 


-7.10 




(T) 


26.69 


39.46 


39.43 


46.68 


46.67 




(Vc) 


-35.06 


-26.26 


-26.24 


-22.00 


-21.98 




(Vboul) 


0.67 


0.64 


0.64 


0.65 


0.65 




( v t) 




-20.89 


-20.88 


-30.47 


-30.47 




(V*) 




-0.03 


-0.03 


-1.97 


-1.97 




V(r 2 ) 


1.74 


1.82 


1.82 


1.79 


1.79 




P(0, 1/2) 


100 


92.98 


92.96 


91.42 


91.41 




P(2,3/2) 




6.98 


6.99 


8.51 


8.53 




P(l,l/2) 




0.03 


0.03 


0.04 


0.04 




P(l,3/2) 




0.02 


0.02 


0.02 


0.02 




4 He(0+) 


E 


-29.94 


-25.29 


-25.26 


-25.08 


-25.05 




(T) 


58.08 


86.93 


86.77 


101.59 


101.36 




(V) 


-88.86 


-66.24 


-66.11 


-54.93 


-54.73 




(Vcoul) 


0.83 


0.76 


0.76 


0.77 


0.77 








-46.62 


-46.55 


-67.85 


-67.79 




<Vb> 




-0.13 


-0.12 


-4.65 


-4.66 




vV> 


1.41 


1.51 


1.51 


1.49 


1.49 




P(0,0) 


100 


88.46 


88.50 


85.76 


85.79 




P(2,2) 




11.30 


11.26 


13.87 


13.85 




^(1,1) 




0.25 


0.24 


0.36 


0.36 





are comparable to state-of-the-art calculations available in literatures. 

We note from the comparison between AV8' and G3RS results that the solu- 
tion with AV8' has large kinetic energy corresponding to the strong short-ranged 
repulsion. The tensor contribution is larger than the central force contribution in 
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AV8'. On the contrary, the G3RS potential gives larger attraction in the central 
contribution than in the tensor contribution. Though they give different poten- 
tial contributions, both lead to almost the same total energy and nucleon rms 
radius. These features of the realistic potentials are very consistent with those 
of Table 1 listed for the deuteron. 

Now we discuss how well the GVR can incorporate unnatural parity compo- 
nents. The unnatural parity state is usually unfavored because it has a larger 
kinetic energy than the natural parity state and its diagonal matrix element 
becomes rather high. If the off-diagonal matrix elements between the unnatu- 
ral parity and natural parity states are large enough to compensate the loss 
of the kinetic energy, the mixing of the unnatural parity state becomes im- 
portant. Table 3 lists the decomposition of the energies of 3 H and a-particle 
according to the (L, S) channel. The value for the diagonal channel denotes 
jTM T \H\ty \ls)JMjTM t ) , while the value between the different 
channels is 2X LS X L , s ,{$ \ L S)JMjTM t \ h \^ \v s>)JMjTM t ) ■ Note that, in the case 
of 3 H, the contributions from the two spin channels of 5=1/2 and 3/2 with the 
unnatural parity of L=l are summed together. The probability of the unnatural 
parity state is negligibly small in 3 H, as shown in Table 2. The probability of the 
unnatural parity state P(l, 1) for a-particle increases to 0.25% for G3RS and 
0.36 % for AV8'. With this mixing of the unnatural parity component the tensor 
coupling with the (22) channel becomes important and thus the a-particle gains 
energy by about 0.8 MeV for G3RS and 1.6 MeV for AV8'. 

The Oj" excited state of 4 He, located at 7.19 MeV below the p+p+n+n thresh- 
old, is a very good example to demonstrate the importance of the unnatural 
parity state. The Oj" state consists of two channels, (L, S)=(l, 1) (natural par- 
ity) and (2,2) (unnatural parity). The matrix element of each operator in the 
Hamiltonian is listed in Table 4 for the 07/ state, with the column-row index 
of the matrix being labeled by the channel, (1, 1) or (2, 2). The values in paren- 
theses stand for the matrix elements, { 1 &(ls)JMjTM t \H\ x S>(l>s').JMjTM t ) etc The 
table confirms that the 0j~ state cannot be predicted well below the p+p+n+n 
threshold if the coupling of the two channels is neglected. The diagonal ma- 
Table 3. Total energies, in MeV, of 3 H and a-particle decomposed into different (L, S) contri- 
butions. 







G3RS 






AV8' 




3 H(f) 
















(0,*) 


(2,1) 


(1,2+D 




(2,|) 


(1,2+D 


(0,*) 


4.09 


-22.54 


-0.00 


9.72 


-33.60 


-0.03 


(2,|) 




10.85 


-0.23 




16.35 


-0.42 


(1,3+1) 






0.11 






0.22 


4 He(0+) 




(0,0) 


(2,2) 


(1,1) 


(0,0) 


(2,2) 


(1,1) 


(0,0) 


0.95 


-46.65 


-0.01 


12.94 


-68.67 


-0.21 


(2,2) 




21.24 


-1.56 




32.31 


-2.90 


(1,1) 






0.72 






1.47 
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Table 4. The Hamiltonian matrix elements, given in MeV, for the 0^ state of 4 He. The column- 
row of the matrix is labeled by the channel (L, S), which is arranged in the order of (1, 1) and 
(2,2). See text for the matrix elements in parentheses. The P(L, S) values are P(l, 1)=0.955 
and P(2,2)=0.045. The G3RS potential is used. 



H 




T 




0.64 


-13.50 


41.19 (43.13) 






6.46 




7.19 (159.7) 


v c 




^Coul 




-27.82 (-29.13) 




0.46 (0.48) 






-1.10 (-24.42) 




0.02 (0.45) 


v t 




Vb 




-13.48 (-14.11) 


-13.51 (-65.15) 


0.28 (0.30) 


0.01 (0.04) 




0.35 (7.81) 




-0.00 (-0.00) 



trix element of the Hamiltonian is 0.64 MeV in (1,1) channel and 6.46 MeV in 
(2, 2) channel. The unnatural parity channel gives large positive energy because 
its kinetic energy is very large. However, the coupling matrix element between 
(1,1) and (2,2) channels amounts to —13.50 MeV, which locates the 0J~ state 
at —6.40 MeV in good agreement with experiment. This large coupling matrix 
element is brought about by the tensor force and its contribution is quite sig- 
nificant. The admixture of the unnatural parity components is so large as 4.5%. 
We repeated a calculation by omitting the (2, 2) channel. The resulting energy 
is approximately — lMeV, which is higher by about 5 MeV than the full chan- 
nel calculation. Though the importance of the tensor force in the 0^~ state was 
pointed out many years ago [27, 28], the present result indicates that not only 
the tensor force but also the kinetic energy and the central force are important 
factors to determine the energy of the 0j~ state. Our calculation ignores another 
unnatural parity component (0,0). However, its contribution is probably negli- 
gible because it has no tensor coupling with the main channel (1, 1). 

It is interesting to ask the following question: How different are the G3RS 
and AV8' potentials? As seen in Table 2, both potentials give similar energies for 
the few-nucleon systems. Figure 1 displays the energy change of the ground and 
first excited states of 4 He as a function of the basis dimension. Here the basis 
set is the one selected using the G3RS Hamiltonian. The basis is increased to 
optimize the 0^~ state after the dimension of 400. In the curves labeled G3RS the 
basis functions are ordered according to the importance criterion for the ground 
state. Now the curves labeled AV8' stand for the energy change obtained by just 
diagonalizing the AV8' Hamiltonian in the same basis set. The energy obtained 
in this way loses only 210 keV compared to the energy listed in Table 2. This 
means that the basis set determined with the G3RS potential already provides a 
fairly good basis set for the AV8' potential. Thus one may skip the basis search 
for AV8' but only needs to fine-tune the basis to reach the converged solution. 
This is particularly helpful in saving computer time because the G3RS potential 
is mild and simple enough to render numerical calculations fast. 

Using Eq. (B.26) and noting the remark below it enables one to calculate the 
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Figure 1. The energy convergence of the ground and first excited states of 4 He with the G3RS 
and AV8' potentials. The basis set is selected for the G3RS potential by the SVM, and the 
same set is used to calculate the energy with the AV8' potential. 



correlation function 



C(r) 



Airr 2 



(5.5) 



which has the normalization of 4ir dr r 2 C{r)=\. The left panel of Fig. 2 plots 
the correlation functions of the a-particle calculated with the different potentials. 
The curve with AV8' agrees very well with that of the benchmark calculation [14]. 
The curves calculated with the realistic forces show a strong depression at short 
distances due to the short-ranged repulsion, which is in sharp contrast with 
the MN potential case. One can see that the AV8' correlation function is more 
strongly suppressed than that of the G3RS potential. All the correlation functions 
are similar at large distances. The right panel of Fig. 2 compares the correlation 



0.025 



0.015 



q- 0.01 



0.005 




.-■•J 


I I I I 

AV8' 




G3RS . 


II \ 

II \ 
ll \ 
It \ 
— ll \ 
It \ 


MN 




i i^-— — j_ — i 



3 

r[fm] 




6 



Figure 2. Left: the correlation functions of a-particle with the different potentials. Right: the 
correlation functions for the ground states of few-nucleon systems calculated with the AV8' 
potential. 
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functions for few-nucleon systems calculated with the AV8' potential. We note 
that the peak position of C(r) is almost independent of the systems. One can 
see that the correlation function of the a-particle has larger amplitude around 
r=lfm and a shorter tail than those of deuteron and 3 H. 

A two-nucleon momentum distribution also gives useful information on the 
correlation. The momentum distribution is easily calculated by transforming the 
coordinate-space wave function to momentum space (see Appendices A and B.3): 

C(k) = (*| 8{\ \k! - k 2 \ - k){P 1N P 2N } + |*) , (5.6) 
where the projection operator of a pair of nucleons, {PinP2n}+, stands for 

PiA \ ( pp\ 

PiuP 2 u for mi . (5.7) 

^(AA + PiA) / \ pn J 

The momentum distribution defined above has the following normalization 

r°° 2 

4vr ^ dk k 2 C(k) = A{A _ x) <*| "£{P iN P jN } + |*) . (5.8) 

We plot in Fig. 3 the np and pp momentum distributions of the a-particle. The 
realistic forces produce similar distributions. The distribution calculated with 
the MN potential follows the realistic distribution up to k < 1.3 fm -1 , but then 
decreases rapidly as the momentum increases. This is because the MN poten- 
tial has neither strong short-ranged repulsion nor tensor force [25]. As shown 
in ref. [29], both np and pp distributions do not show very different behavior 
except that the former is roughly four times larger than the latter as expected 
from the relation (5.8). This is in sharp contrast to a specific momentum distri- 
bution such that the sum of the two momenta, ki + k 2 , is close to zero [25, 29], 
where the nn {pp) momentum distribution, differently from the np distribution, 
is characterized by a dip around 2fm~ 1 . 




Figure 3. Two-nucleon relative momentum distributions of a-particle with the different po- 
tentials. 



Y. Suzuki, W. Horiuchi, M. Orabi, and K. Arai 



15 



Calculations for other J 71 " states of 4 He are in progress. A detailed analysis 
will be reported elsewhere [30] . 

5.2 3a system with a nonlocal RGM potential 

A second test of the CG basis calculation is a 3a system for 12 C. That the act 
potential cannot be local for describing 12 C is not surprising because the in- 
teraction between composite particles is intrinsically nonlocal. Such interactions 
between two composite particles can be derived microscopically within the RGM. 
A first attempt to use an energy-dependent aa RGM kernel has been developed 
in ref. [31] in a study of 12 C. Using some average energy for this kernel provides 
fair results for 12 C [32] but not for other three-cluster systems [33]. Hence the 
energy dependence is eliminated from the RGM equation, which enables us to 
obtain energy-independent nonlocal kernels. Recently, we have shown that this 
procedure leads to a consistent description of the two + states of 12 C [34] as 
well as the 3/2~ and 5/2~ states of 9 Be and the ground state of 6 He [35]. 
The RGM equation for the two-cluster relative motion function \ reads 

(T + V + sK) X = e X , (5.9) 

with a bare RGM potential V 

V = V D + V EX = V D + K T + K V . (5.10) 

Here T is the intercluster kinetic energy, K is the overlap kernel, Vd is local and 
called the direct potential, and the nonlocal potential V EX comprises Kt and Ky 
which are the exchange nonlocal kernels for the kinetic and potential (including 
the Coulomb term) energies, respectively. Equation (5.9) can be converted to an 
energy-independent equation for g = ^/M % = y/1 — K\ as follows: 

(T + V RGM )g = eg (5.11) 

with 

yRGM = ^-1/2( T + V ^-l/2 _ T = V + W ^ ( 5 _ 12 ) 

where W is a difference between the renormalized RGM potential y RGM and the 
bare RGM potential V: 

W = N- l/2 {T + V)N~ 1 ' 2 - (T + V). (5.13) 

The energy-independent potential y RGM serves as an interaction between the 
composite particles. As shown in ref. [35] the nonlocal interaction W can very 
well be approximated with \W^\ where 

W {1) = ^[K(T + V) + (T + V)K], (5.14) 

and A is an appropriate constant. The A value is set to be 1.30 for the aa case. 
The aa potential we use in the present study takes the form V+AW^ 1 ). All of the 
needed aa RGM kernels are explicitly given in ref. [33]. We use the harmonic- 
oscillator parameter b = 1.36 fm for the single-nucleon orbit in the a cluster, 
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and the MN potential [22] with u = 0.94687 as the two-nucleon potential. The 
nonlocal aa potential constructed in this way reproduces the aa phase shifts 
very well. 

To obtain the ground and first excited + states of 12 C, we must have such 
a solution that does not contain Pauli- forbidden states <~p n im- These states are 
defined as redundant states of the 2a RGM equation (5.9). They are 0s, Is and 
Od harmonic-oscillator states. We define an operator r for the 3a system by 

3 

r= r ^ ( 5 - 15 ) 

l>k=l 

where r k \ is a separable nonlocal operator acting on the relative motion of the 
particles k and I: 

r H = r kl (0s) + r kl (is) + r kl {od) (5.16) 

with 

i 

r ki(n£) = I ( Pmm(r k - ri)){ip n i> m {r k - n) | . (5.17) 

m=-e 

Eliminating the forbidden states is performed by adding a pseudopotential 
to the 3a Hamiltonian and finding such a solution that is stable for sufficiently 
large 7 [36]. In practice, the 7 value is set around 10 4 MeV in the present cal- 
culation. Since it is hard to eliminate the forbidden states to high accuracy in 
a nonorthogonal basis set like the CG basis, the use of larger 7 values leads to 
numerically unstable energies. The matrix elements of the nonlocal potentials, 
V, and r{nl), can all be calculated using the formulas in Appendices C 

and D. 

Table 5 compares the present result with that obtained by the HH ap- 
proach [35]. In the HH expansion method the size of the three-body basis depends 
on the maximum hypermomentum K mSuX . The HH calculation in the table uses 
Km&x = 36, and for this value there are 100 HH in the variational expansion of 
the wave function [37]. Also the hyperradius variable is discretized over 35 points 
and thus the Hamiltonian matrix has a size of 3500x3500. In contrast to the HH 
expansion, the solution in the CG approach employs 200 basis functions. The 
energies of both the Of and 0^" states are in good agreement between the two 
calculations. These values agree reasonably well with the experimental energies. 
In particular the energy of the 0^~ state is in excellent agreement with the ob- 
served state called the Hoyle state. The matter rms radius of 12 C is defined by 
adding the empirical value of the a-particle (yvU = 1-479 fm). The radius for 
the ground state agrees between the two calculations. The comparison of the 
HH and CG results indicates that the CG basis is quite flexible to describe the 
different shapes of the wave functions in a small dimension. 

The expectation value of each term of the Hamiltonian is also listed in Ta- 
ble 5. It agrees within 30keV for the Of state between the two methods, but the 
difference is still larger for the 0^~ state. This may be due to that both methods 
obtain the latter state in a bound-state approximation without imposing a proper 
asymptotic condition, though it is located above the 3a threshold. To judge how 
well the elimination of the Pauli-forbidden components is achieved with the 200 
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Table 5. The energy E (MeV) from the 3a threshold and the rms radius sj (r 2 ) (fm) for the 
ground state and the first excited + state of 12 C. The energy- independent aa RGM potential 
is used. The expectation value of each piece of the Hamiltonian is also shown. Experimental 
energies are —7.27 and 0.38 MeV for the 0^ and 0% states, respectively. The inputs are taken as 
h 2 /m N — 41.472 MeVfm 2 , m a —4m N , e 2 =1.44 MeV fm~ 2 , where m a is the mass of a-particle. 



Method 




E 


(T) 


<Vb> 


(V EX + \W^) 




vV> 


CG 


of 


-9.83 


79.06 


-80.34 


-8.55 


7.75 


2.19 




0+ 


0.42 


20.34 


-18.32 


-1.60 


0.68 




HH 


of 


-9.83 


79.09 


-80.36 


-8.56 




2.19 




0+ 


0.43 


21.5 


-19.38 


-1.69 







basis functions, we show the expectation value of T: It is 1.34xl0~ 6 for Of and 
2.19xl0 -7 for O2", respectively. The contribution from the pseudopotential 
is subtracted in the energy listed in the table. 

As discussed in refs. [36, 34], a local aa potential leads to unphysical re- 
sults for the 3a energies even though the potential fits the aa phase shifts. For 
example, the well-known aa potential [38] called the BFW potential gives an 
extremely deep energy of -19.08 MeV for the Of state. This is the mam reason 
why we advocate the importance of using the energy-independent RGM kernel 
as the interaction between the composite particles. Though the BFW potential 
cannot produce a physically acceptable result for the Of state, it may be worth- 
while testing if this potential could be used for setting up a suitable basis set for 
the 3a Hamiltonian with the nonlocal potential. This is so because the BFW po- 
tential whose nuclear part is given by a single Gaussian is by far economical than 
the nonlocal RGM potential from the point of view of computer time. Figure 4 
displays the energy convergence as the basis dimension increases. In this figure 
the basis set is selected to optimize the Of and 0^~ states for the 3a Hamiltonian 
employing the BFW aa potential. This basis set is now just used to diagonalizc 
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Figure 4. The energy convergence of the + states of 12 C with the BFW and nonlocal RGM 
(NL) potentials. The basis set is determined by the SVM for the BFW potential. 
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the Hamiltonian with the nonlocal aa potential. The curves labeled NL display 
the energy convergence of the latter calculation for the two + states. The re- 
sulting energies are — 9.75MeV for the 0+ state and 0.43 MeV for the 0^" state, 
which are close to those of Table 5. The rms radius of the 0^ state is reproduced 
as well. In fact the overlap between the wave functions obtained in the set of 
Table 5 and the common set noted above is 0.9995 for 0^ and 0.9987 for 0+. 
This finding is quite appealing for studying a multi-a system interacting via the 
nonlocal aa potential. We just use the BFW potential to set up a suitable basis 
set and then use the refinement process to adapt the basis set to the nonlocal 
potential. This will be much more economical than selecting from the beginning 
the basis for the Hamiltonian with the nonlocal potential. A calculation along 
this line is in progress for a 4a description of 16 O. 

6 Summary 

We have extended our previous study on the global vector representation to 
cope with the orbital motion with an unnatural parity. The basic idea is to 
introduce two global vectors which are defined by a linear combination of the 
particle coordinates. Combining these angular parts with a spherical part of the 
orbital function provides a flexible basis function. We have shown that the matrix 
elements for most of physical operators of interest including nonlocal operators 
can be derived analytically if the spherical part is taken as a correlated Gaussian. 
The fact that its Fourier transform again takes the same functional form opens 
a wide applicability for studying correlations of a system in momentum space. 

Our numerical test examples included A=3, 4 nuclei interacting via realistic 
interactions of AV8 type. They have demonstrated the applicability and accuracy 
of the present method, in comparison with the partial-wave expansion results. 
We have shown that the 0J" state of 4 He at the excitation energy of 21.01 MeV 
offers a good example to indicate the importance of the tensor coupling between 
the natural parity and unnatural parity components. The correlated Gaussian 
basis was also tested for a nonlocal potential in the 3a system, and it succeeded 
to reproduce the results of the extensive hyper spherical harmonics calculation. 

A correlated Gaussian with the global vectors can easily be adapted to bound- 
state problems of a larger system as the tedious angular momentum couplings 
of partial-wave expansion are avoided. It will be interesting to devise a method 
of applying this flexibility to a continuum calculation of scattering and reactions 
of few-body systems. 

Appendix A: Momentum representation of correlated Gaussian basis 

The Fourier transform of the CG defines the corresponding basis function in momentum space. 
As is shown in Eq. (A. 7) below, it is again a linear combination of the momentum space CG. 
The case with a single GV has recently been carried out in ref. [39]. The momentum space 
CG is useful to evaluate matrix elements such as semi-relativistic kinetic energy or momentum 
distribution between the particles. See Appendix B.3. 

The transformation from the coordinate to momentum space is achieved by a function 



$(fc,cc) = 



1 



ikx 



(A.l) 



(2*0 
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where k is an (N — l)-dimensional column vector whose ith element fej, multiplied by h, is a 
momentum conjugate to a;;. The Fourier transform of the generating function g reads 

(*(fe,a:)|g(AieiUi+A2e2M2;^,a;)) = L_ e i iA ~\ (A.2) 

(det A) 2 

where v = Aieiiti+A 2 e 2 u 2 — ik. Simplifying vA~ 1 v by omitting both Ai and A 2 terms, we have 
— > AiA 2 pei • e 2 — iAiei ■ A~ 1 uik — i\ 2 e 2 ■ A~ 1 u 2 k — ^kA~ 1 k (A. 3) 

with 

p = u^A~ 1 u 2 . (A.4) 

Here the symbol — > indicates that the \f terms can safely be dropped. Substituting Eqs. (A.2) 
and (A. 3) into Eq. (2.9) enables us to derive 

{^(k,x)\F iLlL2)LM (u 1 ,u 2 ,A,x)) = B Ll B L2 - — — -exp (-ifcA^fc 

(det^4) 2 \ z 

min(L^,L 2 ) Ll + L2 _ 2t 

x E ei{Li-e)KL2-ey. p J J de ^ de ^ Y ^ e ^ Y ^ e ^ M 

x(ei • e 2 )\ei ■ l Z %ifc) il_<! (e 2 • A^Hi-ik)^ 1 . (A.5) 

As noted below Eq. (2.10), the above integration has a non- vanishing contribution provided 
the vectors e; couple to their maximum value Li . This can be done using the relation 

l\ 



(ci ■ e 2 r =► ■ B -(-irV2l+T[Y«(ei)^(e2)]oo, (A.6) 

where the symbol =>■ indicates that the angular momentum coupling must be made to its 
maximum value for each e^. The terms (ei ■ A~ 1 uik) Ll ~ e and (e 2 • A- 1 u 2 k) L2 ~ i are expanded 
in a similar manner. Combining these results leads to 



{^(k,x)\F {LlL2)LM (u 1 ,u 2 ,A,x)} 

■Y / ^( L i L '2L-J)(-p) e F (Ll _ eL2 _ e)LM (A- 1 u 1 ,A- 1 U2,A-\k), (A.7) 



(-i) Ll+L2 — 



(detA)l e=0 

where £m = min(Li, L2, \{L\ + L 2 — L)/2]). The coefficient K.(LiL 2 L;£) is given by 



K,(L\L 2 L\ £) = M-r/ +L - Ll - i2 (2£i + l)H(2L 2 + l)!! /(2L 1 + 1)(2L 2 + 1 



(2^+l)!!(2Li-2^+l)!!(2L 2 -2£+l)!! y (2£+l)(2L+l) 
x C{ZL 1 -l;L 1 )C{{.L 2 -l;L2)U{L 1 -lL 1 L 2 ~lL 2 ;tL), (A.8) 

where C is a coefficient which couples two spherical harmonics with the same argument 



C(W 2 ; is) = {2 ^(lfl\] 1) {ll0h0m ' (A ' 9) 

and [/ is a unitary Racah coefficient [9]. 

The Fourier transforms of the CG of Eqs. (2.4) and (2.6) are particularly simple: 

(— i) L 

{$>(k,x)\F {L0)LM (ui,u 2 , A,x)) = 3-F( L0)LM {A~ 1 u i,A~ 1 u 2 , A" 1 , fe), 

(det^) 2 

(®(k,x)\F (L1)LM ( Ul ,U2, A,x)) = ± — '- — -F {L1)LM (A 1 u 1 ,A 1 u 2 ,A \k). (A.10) 

(det.4)2 
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Appendix B: Matrix elements of local operators 
B.l Overlap 

The overlap matrix element between the generating functions becomes 

' (27r )^-i\i 



(fl(A3e3M3+A 4 e 4 «4; A', x)\g(\ 1 e 1 u 1 +\ 2 e 2 U2; A,x)) 



detB 



exp ( ^uB 1 v 



(2tt 



.N-l \ 2 



detS 



j>i=l 



where 



cxp ( ^2 pijXiXjd ■ e 3 j , 

4 

V = ^2 ^i^iUi, 



and 



B = A + A', pij = UiB 1 Uj. 
Following the procedure explained in Section (3) and using Eq. (A. 6), we have 



(B.l) 

(B.2) 
(B.3) 



n 

j>i=l 



(pijei ■ e.j) Vi - 



Pij 



\j>i=l Fl] / L 

X V(pi2Pl2j534P34 L X -pi2 L 2 -pvi L 3 ~p 34 L 4 -p 34 0LL0; LL) 



x [ [Y Ll ( ei )Y L2 (e 2 )] L [Y i3 (e 3 )n 4 (e 4 )] J ^ ■ 
Here the non-negative integer powers pij are restricted by Li as follows: 

Pl2 + PVA + Pl4 = Li , P12 + P23 + P24 = L 2 , 
Pia + P23 + P34 = L 3 , P14+ P24 + P34 = L 4 ■ 

The coefficient X needed for a stretched coupling for each a, 

[y„(ei)y»(e 3 )] 00 [Y b ( ei )Y b (e 4 )] QQ [Y c (e 2 )Y c (e 3 )] QQ [y d (e 2 )y d (e 4 )] QQ 
=> ]T X(abcd- L) [[y a+ 6(ei)y c+ d(e2)] L [Y a+c (e 3 )Y b+d (e 4 )] J 



(B.4) 



(B.5) 



(B.6) 



is given by 



X(abcd;L) 



2L + 1 



(2o+l)(26+l)(2c+l)(2d+l) 
X C(ac;a + c)C(bd;b+d) 



C(ab;a + b)C(cd;c+d) 



a b a + b 
c d c+d 
a+c b+d L 



(B.7) 



where the square bracket [• • • ] stands for a unitary 9-j coefficient [9] , while the coefficient Y 
appearing for the coupling 



[[[y o (ci)n( e 2)],[y c ( e3 )y ( i(e4)]J [[y Q (ei)^( C2 )] A [y 7 (e3)y 4 (e4)] 1 

=► Y(abcdafa6t\vK; LL') [ [y a+ a(ei)y +/3 (e 2 )] L [Y c+1 {e 3 )Y d+s {e 4 )} £ ,] ^ , (B.8) 



is given by 

Y(abcdaf3'y6£\vK; LL') — 



' £ 


£ 


" 




A 


V 


K 




L 


L' 


K 





a b £ c d £ 

a P A 7 S v 

a+a b+p L J [ c+7 d+5 L' 

x C(aa; a + a) C(6/3; b+(3) C(c 7 ; c+7) C(d<5; d+<5). 



(B.9) 
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The overlap matrix element reads 

{F(l 3 l 4 ) L m(U3,U4,A',x)\F( LiL2 ) LM (ui,U2, A,x)) 
X X(p 13 p!4P23P24; L)Y(pi2Pl2P34Pi4 L\ — p 12 L 2 ~Pl2 L3—P34 L4— P34 OLLO; LL). 

As seen from Eq. (B.5), only two of p t j, e.g., p 12 and pi 3 are independent for given Li values. 
The sum in Eq. (B.10) consists of only few terms: just one term with pis = L for the natural 
parity, and three terms with (i) pi 3 =L — 1, p\4 — P23 = 1, (ii) P13 =L — 1, pvi = P34 = 1, (hi) 
Pi3 = L, JJ24 = 1 for the unnatural parity. Other pij values all vanish. 

B.2 Kinetic energy and mean square distance 

The kinetic energy with the center of mass kinetic energy being subtracted is expressed as 
T— T cm = (1/2)7? A-7T, where nj is —ihd/dxj and A is an appropriate (N— l)x (N— 1) symmetric 
matrix. The matrix element between the generating functions is given by 

(p(A 3 e3U3+A4e 4 W4; A', x)\^n A7r|g(Aieiui+A 2 e2«2; A, x)) 



where 



with 



-¥'--M^ L y-'(t" r "')- <bii) 



z = ^2x i e i A i B- 1 u i , R = 3Tr(B~ 1 A'AA), (B.12) 



Ai = A 2 = A', A 3 = A 4 = -A. (B.13) 
The desired matrix element is expressed using the overlap matrix element as follows: 

{F(l 3 l 4 )lm{u3, i*4, A', x)\T - T cul \F (LlL2)LM (u 1 ,u 2 ,A,x)) 

h 2 ( 4 d \ 

= 2 Si 'a — {F(l 3 l 4 )lm(u 3 ,U4, A' ,x)\F (LlL2)LM (ui,U2, A,x)), (B.14) 

2 V j>i=i 0pii ) 

where 

Sn = UiB~ 1 AiAAjB~ 1 Uj. (B.15) 

Since pij appears only as (pij) Pij in (F'\F), its differentiation is elementary. 

A mean square distance is a quantity to characterize the size of a system. What is meant 
by it is the expectation value of an operator such as (r fe — r;) 2 or 5ZfeLi( r fe — &n) 2 ■ Clearly such 
operators can be expressed as xQx with an appropriate symmetric matrix Q. For example, 
Eq. (4.1) enables us to obtain Q = ww for the case of (rt,— r;) 2 . The matrix element of xQx 
between the generating functions takes the form 

(<jr(A3e3M3+A4e4M 4 ; A' , x)\xQx\g(XieiUi+X 2 e2U 2 ; A, x)} 
= (TZ + vB-'QB-'v) (^P) 2 cxp ^5B _1 t)j (B.16) 
with lZ = 3Tr(B~ 1 Q). Comparing to the kinetic energy matrix element (B.ll) leads to 



{F(l 3 l 4 ) L m(u3,U4, A' ,x)\xQx\F (LlL2)LM (u 1 ,u 2 , A,x)) 

d 
dpi. 



= {lZ + 2 ^2 Sij-^— \{F( L3Li)LM {u3,U4,A',x)\F( LlL2 ) LM (ui,U2,A,x)) (B.17) 



j>i=l 

with Sij —UiB~ 1 QB^ 1 Uj. 
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B.3 5-function or multipole operator 

A local potential V(rk—ri) can be expressed with use of rk—ri=wx (see Eq. (4.1)) as 

V(r k -ri) = J drV(r)S(wx-r). (B.18) 

The matrix element of V can thus be obtained through that of the S- function. The single-particle 
coordinate from the center of mass, v^ — xn, is also expressed in terms of Xi with appropriate 
coefficients Wi, and a single-particle operator Z)(r/ S — xm) is represented by the 5-function as 

D(r k -x N ) = j drD(r)8(wx~r), (B.19) 

so that its matrix element again reduces to that of the 5-function. 

The matrix element of the 5-function between the generating functions is given by 

(p(A3e3M3+A4e 4 u 4 ; A',x)\8(wx — r)|#(Aieiwi + A 2 e2W2; A, x)) 

= { detB ) eX V{2 VB v--c(r-wB v) 

(2tv) n - 2 cY ( 1 2 A , , * A \ mnn . 

detB I expl--cr + 2^ PijXiXjei-ej+c^li^iei-rj, (B.20) 

e ' V j>i=l i=l / 

where 

c = ({SB^ 1 wy 1 , 7i = u i B~ 1 w = wB^ 1 u i , p tj = pij — cjijj = UiJ^ 1 ^. (B.21) 
The exponential part including r can be expanded in power series: 

n^^n '-'^ iviww)]. 



i=l 

4 



fl ( 1 )^) 9, T .« j £ v^TT C^iga ; /i) C(g 3 g 4 ; //) C(w*'; t) 
i=i qi J [i/i'e 

x [[[n i (e 1 )y g2 (e 2 )] M [y M (e 3 )F 94 (e4)]^]^^)] 00 - (B.22) 

Using Eqs. (B.4), (B.8) and (B.22) leads to the formula for the desired matrix element 

(F( L3Li ) L 'M'(u3,U4, A', x)\8(wx - r)\F( LlL , 2)LM (u 1 ,u 2 , A, x)) 

( ,\L 1 +L 2 +L+M+L'+M' ( 4 \ , . 



s/W+i \fJi 7 V detB 

(-IgVgM+T,- ( A (-1)« 



py \i>i=l ^ / q, \i=l ^ 

x E {LMlm\L'M')V2l+lX{pviPiiP23P2A\K) 

K\Xfx' lm 

X V(pi2Pl2P34P34 P13+P14 P23 +P24 P13+P23 Pl4 +P24 OkkO; Kk) 

x L 2 -q-2 £3-93 L4-94 q1q2q3q4K.nn' £; LL 1 ) 

x C( qi g2;M)C( 9 3g4;M')C( W '^)^ 1+92+M+94 e'^ r2 y_ m (f). (B.23) 
Here pij and ^ are non-negative integers, which must satisfy the following equations 

P12 +P13 +P14 +qi = Li, P12 +P23 +P24 + 92 = L 2 , 

P13 +P23 +P34 + 93 = L3, Pl4 +P24 +P34 + Q4 = 1/4. (B.24) 

Using the expansion 

S (wx - r) = 15(1 WX J ~ r) E y m (^)y m (f )*, (B.25) 

lm 
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we obtain the matrix element for the multipole operator 

{F(l 3 l 4 ) L 'm'(u3,U4, A', x)\V(\ wx \ )Y tm (wx) \F (Ll L 2 ) LM («i , «2 , A, x)) 



( t\L 1 + L 2 +L+L i / 4 \ /(<) 



(LM^m|L'M'}V27+T rQ- 8 ^ 



c 



V2Z/ + 1 l 1J r * / V detB 

x e ( n '- ir y™^ fa)>« ) e (n SF^) 9 - 

Pii \J>i = l 



3 



x / drr qi+q2+q3+q4+2 e-2 cr V(r) 

x E X (Pi3Pi4P23P24; k) C(<?ig 2 ; C(q 3 q 4 ; y) C(wV) 

Kflfl' 

X y(pi2Pl2P34?334 J513+P14 P23+£>24 J513+P23 P14 +JJ24 OkkO; Kk) 

x Y(Li — qiL2 — q2L3 — q,iL4 — q4qiq2q3q4KH^'i;LL'). (B.26) 

Equation (B.26) with = 00 has many important applications. The matrix element of 
the relative distance \rk — ri\ or rfc — a; at is calculated by choosing V(r) = r with appropriate 
w. The correlation function (5.5) is also easily calculated by choosing V(r)=8(r — a) (with 
a being replaced by r later). Of particular importance in hadron spectroscopy is the matrix 
element for the semi-relativistic kinetic energy, \J (urn) 2 + y? . This matrix element reduces to 
that of the central matrix element with V(r) = \Jr 2 + ji 2 because the Fourier transform of the 
coordinate space CG is expressed in terms of the momentum space CG (see Eq. (A.7)). The 
relative momentum distribution between the particles (5.6) can also easily be calculated in the 
momentum space by choosing V(r) appropriately. 

The spin-orbit matrix element is obtained from the following result 

{g{\3e3U3 + \4e4U4; A' , x)\8(wx — r)(wx x i!;7r)|g(Aieitti + A2e2U2; A, x)} 

x exp {^vB^v-^cir-wB^vf^j , (B.27) 

where z is defined in Eq. (B.12). When the radial form of the spin-orbit potential is scalar, i.e. 
V(r) is a function of r, we may omit c£AB~ 1 wwB~ 1 v thanks to the relation 



/ 



dr V(r)(r x a) exp ( -\c(r -a) 2 ) = 0, (B.28) 



which leads to 

(g(\3e3U 3 + \4e4U4; A' , x)\V (\wx\)(wx x £-n) m \g(\ 1 eiui+ \ 2 e 2 U2; A, x)) 

xexpl--cr 2 + E PijAjA 3 ei • e 3 • + c^7,Aie, • r J , (B.29) 
V j>i=i i=i / 

where (a x b) m stands for — \/2i[a x 6]i m = — (4\/2~7r/3)ia&[Yi (d)Yi(6)]i m . 

Comparing Eq. (B.29) with the matrix element (B.20) for the J-function, it is possible to 
obtain an analytic expression for the spin-orbit matrix element though it is very much involved. 



B.4 Potentials of Gaussian radial form 

When the radial form of V(r) is Gaussian, exp (— c'r /2), the formula turns out to be much 
more concise. Even when its form is not Gaussian, it may be approximated very well by a 
superposition of Gaussians with different c' values. As we show below, the matrix element of 
operators with Gaussian radial form can be expressed using the overlap matrix element. 
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(i) Central potential 
Using Eq. (B.20) yields 

(g(\- A e- A u 3 + \4e4UA\ A', x)\ exp j - ^c'(wx) 2 ^ \g(\ 1 e 1 u 1 + \ 2 e 2 U2; A,x)) 

(^) * - {± h - £M ■ -) • 

Comparing this result with Eq. (B.l) confirms that pij in the overlap is here replaced with 
Pij — (cc' j (c + c'))7i7j. The desired matrix element takes exactly the same form as the overlap 

(F(l 3 l 4 )lm(u 3 , u 4 , A', a;) I exp { - {wx) 2 }\F {LlL . 2)LM {ui, u 2 , A, x)) 

= (—^—7) {F(l 3 L 4 )Lm{U3,U 4 ,A' ,x)\F (LlL2)LM {u l ,U 2 ,A,x))\ , .(B.31) 

The Coulomb potential V(r) = l/r is expressed as a superposition of Gaussian potentials: 

1 = 4= r dte- rH2 . (B.32) 



Thus we obtain 



(3(A3e3«3+A4e4M4; -4', !— r |g(Aieiui + A 2 e2tt2; A, x)) 

\wx\ 

2 r 00 

= —j= I dt(g{\ze A uz + \AeAUA\ A' ,x)\ exp{— t 2 (wx) 2 } |gr(Aienti + A 2 e 2 it2; ^4, a;)} 

/2c f 1 / f27r~) iv_1 \ ' / 4 \ 
— " Y — J du y detB j ex P y Yl [Pv _ c " 2 7i7i] AiAj-e, • ej J , (B.33) 

where a change of the integration variable is performed through t—^/c/2 u/^/l — u 2 . The 
Coulomb matrix element reduces to the following integral of the overlap matrix element 

{F(l 3 l 4 )lm{u3,ua,A',x)\j^^\F^ LiL2)lm (ui,U2,A,x)) 



/2c f 

= \ — / du(F( L3L ^ LM {uz,UA,A',x)\F( LlL2)LM {ui,u2, A,x))\ .(B.34) 

V 7T J() 'Pij^Pij-cu 2 -tifj 

The above integrand is a polynomial function of u whose degree is at most L1+L2+L3+L4, so 
that the integration can be accurately performed with use of the Gauss quadrature. 
We note that the modified Coulomb potential 

I erf t Kr ) = 2lL C d t *r K * rH * (B.35) 
r \px Jo 

can easily be calculated in the same way as above. 

Evaluating the matrix element of Yukawa potential is possible with use of Eq. (B.26). 
Another simple formula is, however, obtained as below. Expressing the Yukawa potential as 

h-"=j*r (b - 36) 

which is a superposition of Gaussian potentials, we get an expression for the matrix element 
(F(l 3 l 4 )lm(u3, «4, A', aO|pJ-r e' KlwXl \F {LlL2)LM (ui ,u 2 ,A, x)) 



^-^.exp^-— j 



x(F i L :i L 4 )LM(u3,UA,A',x)\F (LlL2)LM (u 1 ,U2,A,x))\ . (B.37) 

Pij^Pij-c^lilj 
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As the last example of the Gaussian central potential, we consider V(r) = r exp( — c'r 2 /2), 
which is a derivative of the Gaussian potential, exp(— c'r 2 /2). Though we can evaluate the 
matrix element for this potential from Eq. (B.26), it is possible to relate its matrix element to 
the overlap matrix element. Starting from Eq. (B.20), we obtain 

1 ,,~ 



(sf(A 3 e3'U3 + A4e4'u 4 ; A', x)||ukc| cxp j - -c (wxj 2 ^ |g(Aiei«i + \ 2 e 2 u 2 ; A, x)) 
f (2n) N - 1 Y /A _ . . \ 

) exp y.^ /; A A r ■ 6j j 



tyc \c + c' J \ detB 



( 2 2 4 

X I 1 + [l + ^-^7 lili X i^i e i ■ e j] 

\ 3>i=l 

f 1 f c 2 4 1 \ 

X J dltexpj- — U 2 ) ^2 lilj^i^3 e i ■ e 3 | J ■ (B.38) 



The required matrix thus reads 

1 ,,~ 



(F(l 3 l 4 )lm(u 3 ,U4,A',x)\ \wx\ exp { - -c'(tra) 2 } \F (Li l 2 )lm(ui, u 2 , A, x)) 
= \l \ — 7~ 7 ) f (-P(i 3 L 4 )iM(w3,U4, A',a;)|F( Z , li2 ) i M(iti,U2, A,a;))l 

V 7TC \C + C / y Pij^Pij 

f 1 , r 2c 2 ^ 9 I 

x {F( L - 1 l 4 )lm(u3,u4,A',x)\F {LiL2)lm (u 1 ,u 2 , A,x))\ C(c / +Cll 2) • (B.39) 

Pij ^Pij ^7 lilj J 

(ii) Tensor potential 

For V(r) = exp (— c'r 2 /2) Y2m(r), we express it as Gaussian radial form by using a formula 

e~i c r Y 2m (r) = e-5 c r y 2m (r) / die""- . (B.40) 

Jo 

Then use of Eq. (B.20) leads to 

(5(A3e 3 W3 + A 4 e4«4; A',x)\ exp j - ic'(wcc) 2 jY 2 m(w;c)|.g(Aieiui+A2e2« 2 ; A,x)) 

f(27v) N - 2 cY ( 1 2 x A A A \ 

X ( dctg ) cxp I --cr + 2^ Pi. A. A c. • e 3 • ~ A <' • 

^ ' \ J>i = l i=l / 

x ' i— 1 

X / dt ( c + c' + 2t ) 2 6XP (.£ h + c + C '+2t 7 ^] AlA ^ ' 63 



/ + duu 2 exp ( Ip^ + cw7i7jl AjA^e, • e,- J . (B.41) 



The term y 2m (J2t=i li^i e i) ls expanded as follows: 



4 4 rrz — 4 

^(^T^e,) =^ 7l 2 AX(e ! ) + y / Y 7i7iAiA J -[ri(e i )yi(ej)] 2m . (B.42) 

i— 1 i— 1 j>i=l 
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The exponential part of the integrand of Eq. (B.41) is expanded in power series (AjA^e, • e.j) Vi i 
as in Appendix B.l. Suppose that the pij values give the following angular momenta 

Pl2+ Pl3 + Pl4 = Ll , P12+ P23 + P24 = L 2 , 

Pl3 + P23 + P34 = L 3 , pi4 + p24 + P34 = L 4 ■ (B.43) 

These terms contribute to the matrix element provided Lk — Lk — ^S^i is met in the case of the 
coupling with F 2m (e;), the first term of the right side of Eq. (B.42), while Lk = Lk~8ki~8kj is 
met for the coupling with the second term, [Yi(ei)Yi(ej)] 2m . Let us define the corresponding 
coupling coefficients Z\ and Z2 by 



Y^{Bi) [[Y Ll ( ei )Y- L2 (e 2 )] L [Y- L3 (e,)Y Li (e 4 )] } 
^ ^>i (kLMsLiL, LL'; i) [ [Y Ll ( ei )Y L2 (e 2 )] L [Y Ls (e 3 )F i4 (e 4 )] L ] ^ , (B.44) 



and 



[Fi(e fc )>i(e i )]^[[>L 1 (ei)^ 2 (e2)] z [yz3(e3)y £4 ( e 4)] z ] oo 

=> Z 2 {kU1 2 UUL, LL'; fcj) [ [Y Ll (ei)Yz, a (c a )] L [Yl 3 (e 3 )Y Li (e 4 )] L ,] . (B.45) 

LL' KM 

The Zi coefficients are given by 

z^kUUUul,ll';\) = s LL ,(-iy +L - L I 21+1 



(2k+1)(2L + 1) 



W(kLiL 2 LL), 



Z 1 {kL 1 L 2 L 3 L 4 L,LL';2)=5 Ll 
Z 1 (kL 1 L 2 L 3 L 4 L, LL'; 3) = 



2L + 1 



(2k+1)(2L + 1) 



W(kL 2 LiLL), 



21/ + 1 



(2«+l)(2L+l) 

Zi(kLiL 2 L 3 L 4 L,LL';4) = 



VK(kL 3 L 4 LL'), 



2L' + 1 



(2/t+l)(2L + l) 



with 



W(Kabcd) = C{na; k + a)U(K,adb; K + ac) 
The Z2 coefficients are given by 
Z 2 (kLiL 2 L 3 L 4 L,LL'; 12) 



W(kL 4 L 3 LL') (B.46) 



(B.47) 



= <W-i) 



2L + 1 



-C(1Li;Li)C(1L 2 ;L 2 



(2k + 1)(2L + 1) 
Z 2 (kL\L 2 L 3 L 4 L,LL'; 13) = W (lLiL 2 LL)iy(lL 3 L 4 LL') 



1 1 K 

ii L 2 L 

Li L 2 L 

11k 

L L 

L L' k 



Z 2 (kLiL 2 L 3 L 4 L, LL'; 14) = (-1) Z_L ' +1 W(lLiL 2 LL) W(1L 4 L 3 LL') 
Z2(kL 1 L 2 L3UL,LL'; 23) = (-1) £ ~ L+1 W (1L 2 LiLL) W{IL 3 L 4 LL') 



Z 2 {kLi_L 2 L 3 LaL,LL'; 24) = (-l) L+i ' W {lL 2 L-i_LL)W (1L 4 L 3 LL' 
Z 2 (kLiL 2 L 3 L 4 L,LL'; 34) 



11k 
L L 
L L' k 

l Ik' 
L L 
L L' « _ 

1 Ik" 
L L 
L L' k 



= 6r 



2L' + 1 



(2k + 1)(2L + 1) 



C(1L 3 ;L 3 )C*(1L 4 ;L 4 ) 



1 1 K 

L 3 L 4 L 
L 3 L 4 L' 



(B.48) 
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The tensor matrix element of the Gaussian radial form reads 

( f (l 3 l 4 )l> m> (u:t,v,4., A' , x)\ cxp j - ic' (wx) 2 1 Y 2m (was) \F (L% l 2 ) LM ("I , u 2 , A, x)) 

— - — — — ^=^= — ——{LM2m\L'M') ^ / + 
y/2L' + l 2 y 

mE^- 1 )^ 2 ( II fr 1 ) E ^ 2Z + ^{lUUUUl, LL'; k) 

I fc = l \i=l Zi / Z 

X {F ( l 3 l^lm(u3, u 4 , A', x)\F { i 1 i 2) - L m{u 1 ,u 2 , A, x))\ 



4 



E 7 fc 7i(-l) Ll+L2 [l[^)^V2L + iz 2 (2L 1 L2L 3 L 4 L,LL';kl) 

l>k=l \i=l L i / L 

><{F(L a L i )LM{uz,u 4 ,A\x)\F {LlL2)L M(ui,U2,A,x))\ \. (B.49) 

It should be noted here that Li = Li — 25^ in the first sum of the curly bracket, whereas 
Li = Li — 8ik — 8ii in the second sum. The integral appearing in Eq. (B.49) can be reduced to 

where f(u) is a polynomial function of tt whose degree is at most (Li +L2 + L3 + L4) /2. This 
integral can therefore be accurately evaluated using the Gauss quadrature. 

We note that the calculation of the tensor matrix element of the following radial form 

{F(l 3 l 4 )lm(u3,U4,A', x)\ cxp j - ^c\wx) 2 ^y2 m {wx)\F( LlL2)LM {ui,U2, A,x)) (B.51) 

is easier. In this case the integration of Eq. (B.40) is not needed, and the above matrix element 
(B.51) is obtained simply by replacing u with c/(c + c') in Eq. (B.49). 

(iii) Spin-orbit potential 

Using Eq. (B.29), we obtain 

(p(A 3 e3'U3+A 4 e4U4; A', x)\ exp |-ic'({Sx) 2 | (wx x ^Tr) rn \g(\ 1 e 1 ui+\2e2u 2 ; A, x)) 
^ (^) 2 (^MB-) 2 £ (7^-7^)^7 (efcXei)m 



4 



:exp Yl 

\j>i=l 



CC 



XiXjei ■ e-j (B.52) 



with 



n^iA.B^m. (B.53) 
The difference in the matrix elements between the overlap (B.l) and the above spin-orbit 
potential is that the latter contains an extra factor —i(e.k x e;) m . 

Following the procedure similar to the one described in Eqs. (B.42), (B.43) and (B.45), we 
obtain the spin-orbit matrix element for the Gaussian radial form as follows: 

{F(l 3 l 4 ) L i M > (u3,U4, A' , x)\cxp i^-^c (wx) 2> ^ (wx x £n) m \F( LlL2 ) LM (ui, U2, A, x)) 

( _- ] \L 1 +L 2 + L+L' ( p) \ f 

-(LMlm\L'M') (W^l ( - 



721/ + 1 ' 1 ' \ V 3 J \c + c> 

x ^ kTl - 7^ fc )M) Ll+£2 (nfr- ) E\/2Z + lZ 2 (lZ 1 L 2 L3l4L,LL';fcO 

l>k = l \i=l L i ) L 

x (F(i,£ 4 )£m(«3,B4 1 A' ,x)\F (ilL2)LSl {u 1 ,U2,A,x))\ , , (B.54) 
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where Li is defined by Li — 5ik—5u depending on the summation labels fc and I. 
(iv) Multipole moment 

We calculate the matrix element for a multipole moment operator. Using Eqs. (B.18) and (B.20) 
with V(r) — yim(r), we obtain 

(p(A3e 3 ii3 + A 4 e4U4; A' ,x)\y em (iux)\g(X 1 e 1 ui+X 2 e 2 u 2 ; A,x)) 

dry em (r) ( ^ J expl--cr 2 + ^ p^XiXja ■ e, + 7;A;e; • r J 

^ ' \ j>i=l i=l / 

= f dctg J ^m(^ 7<Aie< ) 6XP ( /'v A A ,; ' ,; • ( B - 55 ) 

^ ' i=l \j>i=l J 

In the case of £ = 2, the above result is compared to Eq. (B.41). That is, it is simply obtained 
by dropping c/2 and by setting u=l. Therefore the matrix element 

<-F(L 3 -t4)iAf(w3, U4, -4', a;)|>'2m(wa;)|F (ili2)iM (Mi, u 2 , A, x)> 

is obtained from Eq. (B.49) by the same procedure as noted above. 

It is easy to obtain the matrix element for the dipole operator (£=1) in a similar way: 

(F(L 3 L 4 )L>M> («3, «4, A', x)\y im (wx)\F {LlL2)LM (u l ,U 2 , A, x)) 

(_X) i i+ i 2+i+i'+i v ^3 



V2L' + 1 

4 / 1 

X 



-{LMlm\L M ) 



J2lk(-lf 1+L2 ( II |^ ) E ^ 2Z + 1 ZxilLMaUL, LL'; fc) 

fe = l \»=1 / Z 

x(-F , (L 3 L 4 )LM(f3,W4, A',x)|F (£i£2)£ m(ui,w 2 ,^, x)>, (B.56) 
where Zi is defined by Li — Li — Sit depending on the summation label fc. 

B.5 Many-particle correlation function 

A many-particle correlation function is useful to visualize the structure of a system [6, 21, 40, 41, 
26]. The function is defined by the matrix element of the product of n 5-functions (n < N—l), 
e.g., 6(ri— Tj— di)5(r k — r t — d 2 ) ■ ■ ■ = 5(w^x — di)S(w l -' 2 '>x — d 2 ) ■ ■ ■ . The basic matrix element 
is the one between the generating functions: 

n 

(g(A 3 e3M3+A4e4W4; A', x)\ Y\ 5(w( a )x - d a )|p(Aieiwi+A 2 e2W2; A, x)) 

a=l 

! ! 1 expf-uB ^--QCQ], (B.57) 



detB J \ (2?r) 

which is obtained by expressing the 5-functions as Fourier integrals, and where 

(C _1 )c«/3 = w&B- l w w , Q a = w(")B~ 1 v — d a , (a, (3 = l,...,n). (B.58) 
Omitting the A? terms leads to 



1_ 



4 4 



2 »B V-iQCQ — > ^ QijXiXjei-ej+^Xiei-Di-^dCd, (B.59) 

J>i = l i=l 

where d is a one-column matrix consisting of di, d n , and 

gij = - (fcr);.;, r>i = (rcd)i (b.60) 

with 

r ai = io(«)B _1 Ui, (a = l,...,n; i = l,...,4). (B.61) 
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Following the same procedure as in Appendix B.3 (cf. Eq. (B.23)) leads to 

n 

(F(L 3 L 4 )L'M'(U3,U4, A' , x)\ Y[ S^^x - d a )\F( LlL2)LM (ui, u 2 , A, x)> 



a=l 



- (n b,) (^) * (m 1 -p 

Pij \J>«=1 PsJ / H \»=1 9l 

x ^ (LM£m|L'M')v/27+TX(pi 3 pi4P23P24;K) 

KfXfx' im 

X y(pi2Pl2P34P34 P13+P14 P23 +P24 Pl3 +P23 P14+P24OKKO; Kk) 

x y(Li — iji L'2 — q2 Li — q3 1/4 — (74 q\qiq3QAKpni ' t; LL 1 ) 

x [^(DiJ^^^^N^^ilj-™. (B.62) 

The values of pij and §i are restricted by the condition of Eq. (B.24). 

Note that the matrix element of a many-body force is readily evaluated using Eq (B.62). 



Appendix C: Integral transform of the correlated Gaussian basis 

C. 1 Integral transform of the generating function 

An integral transform of the CG is needed to adapt nonlocal RGM kernels into a variational 
calculation. A typical form of the kernel we consider is 



TS- ft \ / 1 / II 1 2 I \ 

K(r , r) — exp I — —p r — -pr — qr ■ r I , 



(C.l) 



where r is one of the relative distance vectors, say, between particles k and /. Thus r denotes 
wx (see Eq. (4.1)). The vector r' is defined in exactly the same way as r. 

In order to calculate the integral transform acted by K on the basis function (2.7), we start 
from its action on the generating function 



[Kg](x) = J drK(r',r)g(s;A,x) 



(C.2) 



The operator K acting on the function g changes it to a new function Kg, which is again a 
function of x because r' denotes wx. 

The first step of obtaining [Kg](x) is to transform the coordinate set x to & set of coordi- 
nates y={y l ,-,y N _ l }: 



(i) 



(C.3) 



The vector w (l) is chosen in such a way that the first coordinate y x reduces to wx = rk — ri. 
The choice of other relative coordinates {y 2 , Vn-i}' namely w' 1 ' (i = 2, . . . , N — 1), is not 
unique. Their choice is, however, subject to the condition that any of {y 2 , ...,y N _ 1 } must be 
independent of y lt that is, it may contain and r; only as a combination of (rnkrk-\ j miri)/( / mk+ 
mi) but not as or r; alone; otherwise the integration over r = r^ — r; in Eq. (C.2) cannot 
be performed as it is meant. The coordinate transformation (C.3) from x to y is expressed as 
y = T~ 1 x by the matrix T _1 : 

/ wW \ 



..(JV-l) 



(C.4) 
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The transformation from y to x, x — Ty, reads as 

t= (c (1) C (2) ---C (iV - 1) ), < (1) =C- (C5) 

Since T~ 1 T = TT~ 1 = 1, we have 

JV-l 

w(i ) C U) = C ( i)w 0) = Sijj Cf 'wf ' = Sij. (C.6) 

fc=i 

The second equation is rewritten in a matrix form as 

JV-l JV-l 

^2 (<*>#) = 1 - Cw, or w {k \( k ) = 1 - w(. (C.7) 

fc=2 fe=2 

The next step is to substitute x = Ty to g and to separate the part depending on y x : 
g(s;A,x) = g(Ts\TAT,y) 

= exp^-iayJ + t.y.-^WyW.y^^U^.y^), (C.8) 

where we introduce short-hand notations for the matrix TAT and the vectors, Ts and y, by 

^=(; 1} g), ?-=( t J,), (a9) 

Here a' 1 - 1 is an (N — 2)-dimensional column vector, A^ an (N— 2) x (N— 2) symmetric matrix, 
and is an (iV— 2)-dimensional column vector whose element is a 3-dimensional vector. More 
explicitly they are given by 

a = (AC = ^A(, A$ = C^)A(^ +1 \ 

t = Cs, tf ) =(^)s, (CIO) 

where i, j = 1,2, . . . , N — 2. Substituting Eqs. (C.8) and (C.l) into Eq. (C.2), integrating over 
r — y 1 and then renaming r' ' = y 1 lead to 

[^ ] (,)=(-^-) l e Xp (^-L_^) 5(5 ;A, !/ ), (C.ll) 

where 

/ n' — \ ( 2_+ \ 

M-^> , m -%^y -{*>-%«»)■ <°- i2 > 

The last step is to express the above function g(s; A, y) in terms of x. This is achieved by 
replacing y with T~ x x (see Eqs. (2.11) and (2.12)): 

g(s;A,y) = g^s-T^ AT~\ x) = g{Q K s;A K ,x), (C.13) 

where the matrices Qk and Ak turn out to be given by 

a + p a+p 

A K = A + (a + p - ^—^-)ww - (w(A + A(w) — A(C,A. (C.14) 



a+p / a+p\ ) a+p 

Substitution of Eqs. (C.10) and (C.13) in Eq. (C.ll) leads to the desired result 

[Kg]{x) = {^T~p) 2 GXP i^Tpf^ 3 ) 9(Qks;A k ,x). (C.15) 
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Using this result (C.15) with s = Aiei«i+A2e2W2 in Eq. (2.9), we can derive 



[KF (LlL2)LM (ui , u 2 , A)](x) -- 



2tt 
a+p 



J2K.(L 1 L 2 L ; e) 



a + p 



X F(L 1 -IL 2 -1)Lm{QkUi,Q K U2,A K , as), 



(C.16) 



where £m =min(Li, L2, [(£1 + L2 — L)/2]) and fC(L\L2L; £) is defined in Eq. (A. 8). The action 
of the integral kernel on the CG with two GV leads to that of the same kinds. In particular we 
obtain a very simple result for the CG of types of Eqs. (2.4) and (2.6): 



[KF (L0)LM (ui , u 2 , A)] (x) 



[KF (L1)LM (mi , u 2 , A)] (x) 



2tt 
a+p 

2tv 
a+p 



F(lo)lm(Qku 1 ,Q k u 2 , A k , x), 



F( L1)lm {Qku 1 ,QkU2,A k , x). (C.17) 



Apparently both [Kg](x) and \KF( LlL2 ) LM (ui, 112, A)](x) should not depend on the choice 
of w' 1 ' (i = 2, . . . , N— 1). This is equivalent to the statement that both Qk and Ak are inde- 
pendent of the choice of mjW (i = 2, . . . , N—l). For this to hold true, it is sufficient to show that 
£ is independent of that choice. In order to prove this, let w'^ (i = 2,... ,N—1) denote a set of 
other choice. They are related to the original set (i = 2, . . . ,N—1) by a non-singular linear 
transformation W as 

(V 2 > • • • w'^- 1 ^ = (V 2 > • • • w^) W. (C.18) 
The coordinate transformation from x to y' (y' i = w'( i ^x) is performed by a matrix T" -1 : 



/ w 



(C.19) 



Using Eqs. (C.4), (C.18) and (C.19) enables one to show T'" 1 = ( q ^? ) leading to 

T' — T 



1 _0 

w- 1 



(C.20) 



This equation indicates that the first column of the matrix T' is the same as that of T, namely 
£ , so that C remains the same against any choice of W . 



C.2 Determination of ( 

The vector £ is determined by the condition w(*)£ = Sn. See Eq. (C.6). As shown below, it is 
possible to determine £ without specifying (i = 2, . . . , N—l). 

Any of {y 2 , y N _{\ is expressed as a linear combination of r; as 

JV 

where may be taken arbitrarily provided that they satisfy the following conditions 

N 

mk \ mi \ o rr< oo^ 

a k = ■ A, a t = ■ A, > a t = 0. (C.22) 

m k + mi m k +mi ^ 

The last equation assures that y t , as one of the relative coordinates, has no dependence on the 
center of mass coordinate. The A value is determined as 
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Substitution of n = J2f =1 (U 1 )ijXj (see Eq. (2.3)) enables one to express the vector (C.21) 



in terms of the coordinate x as 



JV JV 

E^ = ^Z( U ~ la )j X j, 
i=l j = l 



(C.24) 



which indicates that the j-th element of is given by (U~ 1 a)j. The condition to determine 
£ is thus expressed as 



E w ^ = E [v 1 )^ - (u-'h] o = i, E ( u ~ la Ui = °- 

3 = 1 3 = 1 3 = 1 

Using Eqs. (C.22) and (C.23) reduces the second equation in Eq. (C.25) to 

N N-l 

E E (^ _1 )«Cj°* = E (^"'OiOi + {u-\) k a k + {U-\) iai 



(C.25) 



i=l 3=1 

= E 

i^k,i 



{u-\)i- 



m k 



(u-\) k 



mi 



en = 0. 



rn k +mi m k + mi 

Since a t (i =fc k, I) can be taken arbitrarily, the coefhcients of a t must vanish, that is, 



N-l 

E 

3 = 1 



{U-%3 



mi 



Cj = for i ^ k, I. 



(C.26) 



(C.27) 



m k + mi * ' " J m k + mi 
The solution of the first equation in Eq. (C.25) together with the above equation determines £. 

C.3 £ 2 -dependent potential of Gaussian radial form 

As a useful application of the decomposition of Eq. (C.8), we calculate the matrix element of 
an .^-dependent potential of Gaussian radial form, where the angular momentum I is defined 
as i = y 1 x (— ih)d/dy 1 with y 1 =wx = r k — n. Using Eq. (C.8) we obtain 

I g(s; A, x) = ^wx x (t - aWy (1) j g(s; A, x) = -ra x (<fs - C,Ax) g(s; A, x), (C.28) 

where use is made of t = C,s and aWj/ 1 - 1 ~C,Ax — awx in the last step. 

Using the property of 1 2 =£ ■ I = $ ■ I and the fact that I commutes with a scalar function 
leads to an expression for the matrix element between the generating functions as follows: 



(g(s';A',x)\exp^~c'(wx)' 1 }■ ( 2 \g(s: \.x)) 
= h 2 {g{s'-,A',x)\cxp\-^c'(wxY \ C\g(s;A,x)) 



(C.29) 



with 



C 2 = (wx x (fa' - C,A'xfj ■ (wx x (Cs - C,Ax)^j 

— xwwx (s'CC s — s'C,C,Ax — s^A'x + xA' C,C,Ax^ 
— (s'C«5a; — xA'(wx) (s(wx — xAC,wx) . 
The integration in Eq. (C.29) is elementary. It reduces to the following result 

( g(A3e3W3+A4e4«4; A', x) \ exp 1 — ^-c 1 (wx) 2 \ t 2 \g(\ieiu 1 + \2e2U2; A, x)) 



(C.30) 



c + c' 



(2n) 



detB 



c + d 



77,7, 



exp ( E Pa 

\3>i=l - 
4 4 

x ( L (0) + 2 E L^fXiXjei ■ e 3 ■+ 4 E E L l?M XiX i ei ' e J X khe k • c, ) . (C.31) 

3>i=l j>i=ll>k = l 
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Here the constants L^°\ and L^ k[ are determined through w, f, Mi, «2, «3, u±, ^4, -A', c'. 
Their expressions are lengthy and not written here. Substitution of Eq. (C.31) into Eq. (3.2) 
makes it possible to relate the desired matrix element to the overlap matrix element by 

(F(l 3 l 4 )lm(u 3 ,U4, A',x)\ exp j - {wxf}t 2 \F (LlL2)LM (u 1 , u 2 , A, x)) 

x(F( L3L4 ) L M(u s ,U4,A',x)\F iLlL2)LM (u 1 ,u 2 ,A,x))\ , . (C.32) 



Appendix D: Matrix elements of nonlocal operators 

D.l Nonlocal operator matrix element for the generating function 

To prove Eq. (4.4), we start from the matrix element of the nonlocal operator (4.3) for the 
generating function. With use of Eq. (C.8) on the bra-ket generating functions, we obtain [21] 

{g(s; A' , x) | Vnl | g(s; A, x)}= J J dr' drV{r' , r) exp ^~^ a ' r ' 2 + *'' r ' — \ ar2 *' 



x I detT 



(2«) 



det{A<- 



^ Ijy ) 5 exp(l^) + ^) ) - d ), (D.l) 



where |detT| 3 is the Jacobian, det(dxfdy), for the transformation from x to y, and 



d = t™ -a«r + t'« -a''V = 



^3) 

V J 



D, 



with 



(D.2) 



(D.3) 



D = s + s' - A£r - A'Cr'. 

Here a', a'' 1 ^, A'^\ t' , t'' 1 ' are defined in exactly the same way as the corresponding quantities 
given in Eq. (C.10) by replacing A and s with A' and s' . 

The matrix element (D.l) must be independent of the choice of (w^ 2 \ . . . , w/" -1 ') from its 
construction. This is indeed assured by a theorem 



det(A (1) + A' w ) = i (detT) 2 detB, 
d(A w +A' w )- 1 d = DJ- 1 D, 



(D.4) 



where c and J are given in Eq. (4.5). Substituting these results in Eq. (D.l) leads to the 
result (4.4), which is manifestly independent of the choice of kj' 1 '. 



D.2 Proof of a theorem (D.4) 

With a slight change of notation the above theorem is stated as follows. 

(Theorem) Let A be an (N — 1) x (N— 1) real, symmetric, positive-definite matrix and T be an 
(iV — 1) x (N— 1) non-singular matrix given by Eq. (C.5). Let a matrix TAT be decomposed into 



TAT = 



a aW 



(D.5) 
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where a, and are defined by Eq. (C.10). The theorem reads 

detA w = i(detT) 2 detA, (D.6) 

/ cm \ 



^(2)^(3) . 



CCD 



A 1 — aA 1 wwA 1 , 



with 



a = (wA 1 w) 1 , 
and a row vector w is the first row of T^ 1 as defined in Eq. (C.4). 



(D.7) 



(D.8) 



(Proof) The heart of this theorem lies in that the quantities on the left sides of Eqs. (D.6) and 
(D.7) are determined by w alone. The theorem was proved by one of the authors (Y.S.) [21]. We 
here show its proof as the application of nonlocal operators becomes increasingly important. 
We can make use of the following identity 



a aW 



1 -iaM 

a 

1 



a 

a w A m-± a m^X) 



(D.9) 



because a = (A( is not zero thanks to the positive-definiteness of A. Taking the determinant of 
the above equation, we have 



det 



a aW 



= (detT) 2 detA = a det [A w - 
Using the Sherman-Morrison formula [9, 42] reduces the right side of Eq. (D.IO) to 

a(l-i^)( J 4W)- 1 a w )det^ 1) . 

To calculate a; = a( 1 )(A (1) ) _1 a (1) , we invert Eq. (D.9), obtaining 

1 ^ a(1) | T~ 1 A~ 1 f~ 1 
1 



(D.IO) 



(D.ll) 







a \ a J \ a , 

(D.12) 

Let G denote an (N — 2)x (N—2) matrix which is obtained by removing the first row and first 
column from the above matrix. Comparing both sides of Eq. (D.12) leads to 



V 



A 



-1 ( (2) (3) (N-l)\ 
[ w 'w v ' ■ ■ ■ W y ' j 



= ( A W-±aV a m) =(AW)- 1 + —(A m r 1 aWam(AV)-\ (D.13) 



where the Sherman-Morrison formula is used in the last step. Multiplying at 1 ) from the left and 
a' 1 ' from the right, we obtain 



T7\^ m x ax 

aWGa (1> = x H = . 

a — x a — x 

Using the definition of G and Eq. (C.10) for makes it possible to obtain 

. 2 

aWGa {1) = CA(1 - Cw)A _1 (l - wQAC, = a 2 wA~ 1 w - a = — -c 



(D.14) 



(D.15) 
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owing to Eqs. (C.6) and (C.7). Thus we obtain x—a~ a. From Eqs. (D.10) and (D.ll), we have 
det^ (1) =(l/(a-:r))(detT) 2 detyl, which is the first equality (D.6). 

To prove the second equality (D.7), we start from Eq. (D.13), yielding the equation 
G~ 1 =A^ — ia'^oW. Inverting with use of the Sherman-Morrison formula leads to 



(A^y 1 = (G- 1 + -a^aWy 1 = G J, Ga^a^G, 

a a + aWGaW 



(D.16) 



where a+aWGa^ = a 2 /a. The left side of Eq. (D.7) is thus expressed as 



£(3) 



(c (2) c 



( 3 ) . . .^-i) ) g 



\ c (JV - l} / 
_ii( C (2) C (3)... C (^-i)) Ga d)^ G 



( <:jv \ 

C (3) 

V J 



((3) 

V J 



(D.17) 



By using Eqs. (C.6) and (C.7) together with the definition of a and a, the first and second 
terms on the right side of Eq. (D.17) turn out to be 

(1 - CwU _1 (l - w() = A' 1 - (QA- 1 - A^wC, + -(C, 

a 

~{1 - C^A-'il - wC)ACCA(l - tw)A-\l - wC) 
a 2 - 

= -- ^ ( A' 1 - (wA- 1 - A- 1 wC+ -CC I ^CC^4 ( A" 1 - C™^ 1 - A _1 wC + — CCJ 
or \ a / \ a / 

= -a ( -ArCC ~ -^ _1 wC - ^Cw^ 1 + A" 1 ™^" 1 ] . (D.18) 
\cr a a J 

Putting together these results completes the second equality (D.7). 
D.3 Examples of nonlocal kernels 

We here collect the matrix elements of typical nonlocal operators which appear in the RGM 
treatment for light nuclei. See Eq. (4.3) for the definition of the nonlocal operator. 

(i) V(r', r) = exp (— \p'v 12 — \pr 2 — qr' ■ r) 

The nonlocal operator is nothing but K(r',r) of Eq. (C.l). Utilizing Eq. (C.16) we obtain 



{F(L 3 L 4 ) L 'M'{U3,U4,A',x) \ VnL \ F(L 1 L 2 )1M ("1 , 1*2, A, x)) 

= (F ( L 3 L i )L'M'(u3,u 4 ,A',x) | [KF (LlL2)LM (u 1 ,u 2 , A)](x)). 
As another method, it is possible to start from Eq. (4.4), yielding 

(p(A 3 e 3 u 3 + A 4 e 4 it4; A', x)\ VNL|s(AieiMi +A 2 e 2 it2; A, x)) 

v / - — '— exp -v.J 1 v + -ZM 1 Z 
detB ) \detM ) 1 \2 2 

where M is a 2x2 symmetric matrix defined by 

ot + p P + q 



(D.19) 



M = 



13 + q a'+p' 



(D.20) 



(D.21) 
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with a, a' , (5 being given in Eq. (4.6) and Z is a 2-dimensional vector whose elements are usual 
vectors defined by 

4 4 

Zi^C-»fl«-'J«' = X>iei/i (i) , Z 2 = (C-'?)s'-'?s = ^A, e ,/ 2 (l) . (D.22) 

i=l i=l 

Here and /j 1 ' are components of a 2-dimensional vector /W; 

/« = f ^ ~ ^ 1 for i = 1, 2, /« = f ) for i = 3, 4. (D.23) 

y -rftii J V C«i - »? «i / 

The desired matrix element is expressed in terms of the overlap matrix element (B.10) as 

{F(L 3 L 4 )L'M>(U3,U4,A',X) | U N L | F {LlL2)LM {u 1 ,U 2 , A,x)) 
3 

/ 27TC \ 2 

VdetAjJ < i? (£ 3 £4)i'M'(M 3 ,W4,^',a;) j F (ilI , 2 ) iM («i, U2, A, x)>| , (D.24) 

where 

Xa = UiJ^Uj + /COM- 1 / ' = Pij + /WM" 1 / '. (D.25) 
As a straightforward application of the above formula, we consider a nonlocal operator 
V(r', r)={erf (r | cr+cr'r' |)/| ar + u'r' |} cxp (— p'r' 2 /2— pr 2 /2 — qr' ■ r) which often appears 
as the RGM kernel for the Coulomb potential. Using the relation (B.35), the operator is ex- 
pressed as an integral of the nonlocal operator of type (i), so that its matrix element reduces 
to Eq. (D.19) with the replacement of p, p' , q by 

p^> p + 2a 2 r 2 t 2 , p ->p +2a' 2 T 2 t 2 , q -> q + 2a& 'r 2 t 2 . (D.26) 

(ii) V(r', r) = (Qur 2 + Q 22 r' 2 + 2Q 12 r ■ r') cxp {-\p'r 12 - \pr 2 - qr' ■ r) 

Using Eq. (4.4) leads to the following expression 

(g(A 3 e3W3 + A 4 e4U4; A' ,x)\VNh\g(^ieiu 1 + \ 2 e 2 u 2 ; A,x)} 
= ^3TrM _1 Q + ZM~ 1 QM~ 1 Z^j 

x (^f^) 1 (H) 1 -P (Uj-^-ZM-^z) , (D.27) 
V dctB J \detM J 1 \2 2 J' v ; 

with 

a -{& £)' < D ' 28 > 

We can follow the steps of Eqs. (B.l), (B.ll) and (B.14), and find that the desired matrix 
element is given by 

(F(L 3 L i )L'M'(.U3,U4,A',x) | VnL | Ffi^jJLM «2, A, x)> 

"(ot)' + ' t 

x(-F 1 (L 3 i 4 )i'M'(M3,M4, | F( LlL2 ) LM (ui,U2, A,x)}\ , (D.29) 

where 

fe 3 = /WM-'QM" 1 /''. (D.30) 

(iii) U(r', r) = (-l/^/27 + T[;^M*• , ):^Mr^)]oo e*P (-^V 2 - |pr 2 ) 

As the last example, we consider a separable kernel which often appears in eliminating redun- 
dant states from the relative motion between the particles. Expressing [3M r ')3M r )]oo as 

PMr')tt(r)]oo = / de'de[y,(e')^(e)]oo^^exp(Ae.r + A' e '.r')| A=v=o , (D.31) 
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{g(\ a e a u a + \4e4U4; A', x)\VNh\g(^ieiui+ \2e2u2; A, as)) 
= (-l)WTT (§) 2 1 1 de'rfe[F,(e')^(e)]oo^ 

v / - ; ' exp I -vJ v + -WM W)\ 

detB ) \detMoJ 1 \2 2 J \\=\>=o 

,£ 

V7>»=1 



: I J de'de[Y e (e')Ye(e)]oo-^r e 
: exp ( V AAiX ( i°e ■ e, + V A'WjjV • e< + (Afo^^AA'e • e' ) I , (D.32) 



where Mo is obtained by putting <? = in M of Eq. (D.21), and 2-component vectors an d 
W are defined by 



X (<) =M -V (<) (i = l,...,4), W=( £ + j (!).:«) 



with Zi and Z2 being defined in Eq. (D.22). 

The desired matrix element is obtained from Eq. (3.2). To simplify the needed calculation, 
we introduce the following notations 



L5 = L§ = £, e 5 = e, e& — e', Ag = A, Ae = A', 

X«=X ( i\ Xa = X { ?, X56 = (Mo- 1 )^, (D.34) 



which make it possible to extend Eq. (3.2) to 



{F(L 3 L 4 ) L 'M'(U3,U4,A',x) I VnL | ^i^jJLM^l, 1*2, A, x)> 

-(-™(n^)(^) , ( 3 ^)* 



n dei([F i3 (e3)yL 4 (e 4 )]^M') a '[^ 1 (ei)^ 2 (e2)] i M[n 5 (e5)F i6 (e 6 )]oo 
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The above operation can be performed similarly to the overlap case, yielding a formula 

{F(L 3 L 4 ) L 'M'('U S ,U4,A',x) \ Vnl j F (LlL2)LM (u 1 ,U2,A,x)} 



'LL'°MM 



(-1) 



L 1 +L 2 +t 



21+1 



2L + \ 



(27T) 



detS 



(JULY 

V detMo J 



E ( n ^/z^ ixvr 



Pij \j>i=l 



Pl2 


P12 





L\—p\2 


L2—pi2 


L 


Li 


L 2 


L 



(2L+1)(2L 5 + 1) 
(2p 5 6 + l)(2L 5 -2p 56 + l) 

X C(pi2 il-pi2;il)C(pi 3 +pl4 P15+P16! L 1 -p r2 ) 

X C(pi2 i2— P12; L2) C(p23+P24 P25+P26; ^2 — P12) 

X C(p 3 4 i3-P34;i3)C(pi 3 +P23 P35 +P36! L 3 -J?34) 

X C(p 3 4 I/4-P34; 1/4) C(pi4+P24 P45+P46', L4-P34) 

X C(p56 i5-p56; is) C(pi 5 +P25 P35+P45; L 5 -p 5 (i) 

X C(p 5 6 L6-P56; is) C(pi 6 +P26 P36+P46; L6-P56) 

xV - 1 

^ V /(2A+1)(2A' + 1)3 



P34 P34 

1/3— P34 1/4— P34 i 
Z/3 Z/4 Z/ 



^(Pl3Pl4P23P24; A)X(pi 5 pi6P25P26; A')X(p 3 5P36P45P46; A') 



P13 +P23 Pl4 +P24 A 
P35 +P36 P45 +P46 A' 
L3—P34 L4— P34 L 



where non-negative integers p^ satisfy the following equations 

P12 + Pl3 + P14 + P15 + Pl6 = Ll, P12 + P23 + P24 + P25 + P26 = L 2 , 
Pl3 + P23 + P34 + P35 + P36 = L 3 , pi 4 + p 2 4 + P34 + P45 + P46 = L4, 
Pl5 + P25 + P35 + P45 + P56 = L 6 , P16 + P26 + P36 + P46 + P56 = L 6 . 





P13+P14 


P23 +P24 


A 


X 


P15+P16 


P25 +P26 


A' 




L1-P12 


1/2— Pl2 


L 




P15 +P25 


P16+P26 


A' " 


X 


P35 +P45 


P36 +P46 


A' 




L5— P56 


^6— P56 






(D.36) 



(D.37) 



There are 15 pij, which satisfy 6 conditions. Thus we have in general 9 independent p„ in 
contrast to 2 in the overlap case. In the case of natural parity states, however, because of 
L2 — L4 — O only two of p^, e.g., P13 and pis are independent. 
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